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ABSTRACT 

It is shown that so-called dark energy could possible be a manifestation of the gravita- 
tional vortex producing the 'gravitomagnetic' (GM) force field: associated with cosmic 
matter rotation and inertial spacetime frame dragging. The general relativistic Godel- 
Obukhov spacetime metric which incorporates expansion and rotation of the Universe 
is used to evaluate this force. This metric is expressed here in spherical comoving co- 
ordinates. Through a cosmic time-scale evolution, it is shown that cosmic acceleration 
is expected when the magnitude of the radial repulsive GM force exceeds that of the 
usual attractive gravitational 'gravitoelectric' (GE) force: associated with just cosmic 
matter and spacetime warping. In general, this phenomenon of cosmic accelerated ex- 
pansion appears to have occurred twice in the history of the Universe: the inflationary 
phase and the present-day acceleration phase. It is suggested in this model that the 
two phases may or may not be related. The cosmological model presented here is de- 
scribed in the context of Einstein's Theory of General Relativity in Riemann-Cartan 
spacetime, which includes cosmic rotation, its effects, and it being considered as an 
intrinsic part of spacetime. Also, a derived analytical expression for the cosmic pri- 
mordial magnetic field is presented; and how it might relate to the GM field through 
the spin density of the cosmic matter is discussed. 

Key words: cosmology: miscellaneous - cosmology: theory - dark energy - gravita- 
tion ~ methods: analytical - methods: numerical. 



1 INTRODUCTION 

Einstein's Theory of General Relativity together with ordinary matter, described by the standard model of particle physics, 
cannot fully explain the observational data from Type la supernovae (Perlmutter et al. 1998; Perlmutter et al. 1999; Riess 
et al. 1998; Riess et al. 2004), the matter power spectrum of large scale structure (Tegmark et al. 2004), and the anisotropy 
spectrum of the cosmic microwave background radiation (Spergel et al. 2007), with all these data suggesting the presence 
of 'dark energy.' From general relativity, assuming homogeneity and isotropy, the standard cosmological model is commonly 
described by the Friedmann, Lemaitre, Robertson, Walker (FLRW) 1920s and 1930s solutions to the Einstein field equations 
for an expanding universe (equation I32|) . According to this standard cosmological model, the expansion of the Universe, if 
it contains only non-negative mass-energy density p and pressure p, decelerates, as expected on the grounds that gravity 
is attractive and the cosmological constant A is zero. The recent observations of cosmic acceleration, first discovered and 
confirmed by Perlmutter et al. (1998) and Riess et al. (1998) from type la supernovae, can only be explained by considering 
repulsive gravity. In the standard cosmological model, this is achieved by models introducing matter with negative pressure 
and/or A ^ [see Perivolaeopoulos (2006), and references therein]. Observations suggest that the alleged cosmic acceleration 
can only be a very recent phenomenon and must have set in during the late stages of the matter dominated expansion of the 
Universe. Recent observations by Riess et al. (2004) identify the transition from a decelerating to an accelerating universe 
to be at 2 = 0.46 ± 0.13. Now, among suffering from the coincidence problem (e.g., Why is the energy density of matter and 
radiation nearly equal to the dark energy density today?) and the cosmological constant problem (e.g.. How could A have 
been so large during infiation but so incredibly small today?), the standard model does not explain why the acceleration has 
started in the recent past. 
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To avoid resorting to anthropic principle arguments to gain acceptance of the above mentioned models, of p < and/or 
A 7^ 0, which are constrained by the so-called standard cosmological model, perhaps we should seek a wider understanding of 
a general relativistic cosmology, not constrained by non-rotation, a spacetime being an extension of the standard cosmological 
model. This 'new' standard cosmological model then should take into account cosmic rotation as well as cosmic expansion: 
two degrees of freedom. When this is done, we find that a repulsive force of gravity is a natural occurrence and could 
possibly provide an explanation for the recent acceleration phase of the present-day Universe, and possibly shed light on our 
understanding of the physics of so-called inflation in the early Universe. 

In this paper, the nature of so-called dark energy is investigated. The aim is to answer the question. Could dark energy 
be a manifestation of gravity? i.e., Could it be that component of gravity, the so-called gravitomagnetic (CM) force field, 
associated with cosmic rotation and inertial spacetime frame dragging? Inertial is used here in the general dynamical sense. 
The Godel-Obukhov metric (Obukhov 2000; Jain, Modgil, & Ralston 2007) which incorporates expansion and rotation of the 
Universe, derived from general relativity, is used to define spacetime separation (or distance), in this quest to answer the 
above question. Importantly, the Godel-Obukhov metric or geodesic line element can ensure the absence of closed timelike 
curves, making it completely causal, different from the originally proposed Godel metric (Godel 1949). The Godel-Obukhov 
cosmological model (Obukhov 2000) contain parameters which smoothly interpolate between this cosmology and the standard 
FLRW cosmology (which describes an isotropic and homogeneous universe filled with matter: commonly represented by an 
ideal fluid). Note, the independent nature of vorticity as associated with shear of a fluid and pure rotation does not allow limits 
on cosmic rotation to be placed by limits on vorticity (Obukhov 2000; Obukhov, Chrobok, & Scherfner 2002; Jain et al. 2007). 
Namely, the Godel-Obukhov spacetime metric is shear-free but the vorticity and expansion are nontrivial. It is not vorticity 
of pure cosmic rotation that would lead to anisotropy of the microwave background radiation temperature distribution, but 
effects of vorticity associated with a shearing force. 

The Godel-Obukhov model does not conflict with any known cosmological observations. Among distinctive predictions 
of the model are effects on the propagation of light (Jain et al. 2007). Cosmic rotation affects a polarization of radiation, 
which propagates in this curved spacetime, resulting in some observable anisotropy (Obukhov 2000). This appears to have 
been confirmed (Birch 1982, 1983; Kendall & Young 1984; Nodland & Ralston 1997a,b; Jain & Saralo 2006). Observational 
tests have been done, that do not require redshift information, by Jain & Ralston (1999). They found significant signals of 
anisotropy in a large sample of data. Several other observations of radiation propagating on cosmological scales have been 
found to indicate a preferred direction, all of which are aligned along the same axis (e.g., Ralston & Jain 2004; Hutsemekers 
et al. 2005). The origin of these effects, however, may be independent of gravitation and restricted to modifications of the 
electromagnetic sector in which polarization observations are exquisitely sensitive (Jain et al. 2007). Perhaps, the model 
presented in this present manuscript will lend support to the possibility that such effects may indeed be gravity related. 

Dark energy is the popular motivation to consider models beyond the standard Friedmann- Robertson- Walker spacetime 
metric. Current data (Jarosik et al. 2011; Tegmark et al. 2004) appear to fit a fiat cosmology with ilmat ~ 0.27 and Q.a ~ 0.73 
for matter and dark energy density parameters, respectively, in the popular lambda cold dark matter (ACDM) cosmological 
model, which assume negative pressure in a FLRW cosmology. These fits assume an isotropic universe, while, at face value, the 
data used in the fits substantially contradicts isotropy (Jain et al. 2007), at least it appears so from the above observational 
tests. Jain et al. (2007) used large redshift type la Supernova data (see Jain et al. 2007, and references therein) and related 
magnitudes, to place constraints upon parameters appearing in the Godel-Obukhov metric (which does not have the restriction 
of an isotropic universe). This is done by obtaining bounds on an anisotropic redshift versus magnitude relationship and on 
accompanying parameters of the Godel-Obukhov metric. They found that the outcome depends on what are used for the host 
galaxy extinctions. The most reasonable fits do not show any signals requiring anisotropy. Yet, the existence of some small 
anisotropy cannot be ruled out. It appears that their findings are consistent with present-day observations, and it might be 
reasonable to investigate models that perhaps yield some anisotropy, particularly the Godel-Obukhov model. 

The Godel-Obukhov metric, with exact general relativistic solutions as expressed by the cosmic scale factor R — R{t) 
and its derivatives (Obukhov 2000), describing the evolution of R{t), just as commonly done for the Friedmann- Robertson- 
Walker metric, avoids the principal difficulties of old cosmological models with rotation, where R{t) describes the expansion of 
physical spatial distances. For example, the Godel-Obukhov metric is consistent with isotropy of the microwave background 
radiation, like the standard cosmology, and it produces no parallax effects. The final state according to this metric depends 
on the values of two cosmological coupling constants (discussed below) in which torsion can cause the Universe to either 
accelerate (■£ > 0) or prevent cosmological collapse (-^ ~ 0), with these constants playing a similar role to that played by the 
cosmological constant. A, in the FLRW spacetime cosmology. 

Using the Godel-Obukhov metric to define separation of spacetime events, we find that cosmic acceleration is expected 
when the radial repulsive GM force (associated with rotational energy) exceeds the attractive radial 'gravitoelectric' (GE) 
force (associated with rest mass energy or mass-energy) This appears to be somewhat the idea behind Einstein's introduction 

^ The terms gravitomagnetic and gravitoelectric defined in this manuscript are not the same as those defined in so-called gravitoelectro- 
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of the cosmological constant A, when he introduced it in his Theory of General Relativity to explain how the Universe could 
resist collapse under the inward force of gravity. However it was later thrown out by Einstein as his greatest blunder. It seems 
that he did not see any use for A in a universe already expanding according to the Hubble law, before discarding it. In any 
case, maybe A, in a general sense, is like the Hubble constant that changes over time with the age of the Universe. That is, just 
as the Hubble constant relates to cosmic expansion of the Universe, and, as we shall see, cosmic rotation, perhaps A relates 
only to the effects of cosmic rotation. In the Godel-Obukhov spacetime cosmology, presented here, two cosmological coupling 
constants [Ai and A3 (Obukhov 2000)] due to rotation and torsion [as related to general relativistic frame dragging or the 
so-called Lense-Thirring (1918) effect] in curved spacetime take on the role of A. Moreover, we find that the Godel-Obukhov 
spacetime cosmology appears not only to explain recent observations of the accelerated expansion, but suggests a dynamical 
description of the Universe over time that is a natural general relativistic extension of the standard FLRW cosmology. In this 
sense the answer to the question posed in the title appears to be yes. 

The organization of this paper is as follows. In Section [2l a detailed description is presented of the model used here to 
explain the present epoch acceleration of the Universe: as being a gravitational-rotational-inertial phenomenon. A formalism 
containing the astrophysical and mathematical descriptions of the components to validate the model's claims is presented 
in Section [3] It includes analytical derivations of the Godel-Obukhov metric in spherical comoving coordinates, the cosmic 
radial GM and GE force fields, and the density of the Universe, where the GE and GM fields are the gravitational analogues 
of electric and magnetic force fields, respectively. In general, the GE and GM fields relate directly to the total mass and 
rotation, respectively, of a gravitating system (see Williams 2002). Also, included are the cosmological parameters used, which 
includes the cosmic rotational (or angular) velocity, the scale factor, and the Hubble constant. The numerical results from 
evolving the analytical expressions for the GE and GM accelerations over time are presented in Section |4l and a discussion 
is presented in Section [S] Included in Section [SI an analytical expression for the cosmic magnetic field (equation 74) given 
by the Godel-Obukhov metric in terms of the spin density is evolved over time and compared with observations and theory. 
Surprisingly, this expression suggests how the magnetic field and the GM field might be related through the spin density. 
Conclusions are presented in Section [G] 



2 MODEL DESCRIPTION 

It seems reasonable to assume that at t ~ the Universe had, at least, two degrees of freedom: translational, associated with 
the expansion and collapse (or infall) in the er-direction and rotational, associated with cosmic rotation in the e^-direction, 
where and are global unit vectors connected to the cosmological spacetime metric. Rotation at t « is consistent with 
the Godel-Obukhov spacetime metric (Obukhov 2000), which allows for rotation and expansion. 

For the initial conditions of the Universe we will assume that a gravitationally unbound 'hot' energetic rotating plasma 
existed at t ~ 10~*^ s, the Planck time, corresponding to the Planck length, Ip ~ 1.6 x 10~^^ cm, wherein beyond it is believed 
that quantum-gravity must be used, since both general relativity and quantum mechanics simultaneously become important 
on such small scales and high density. We assume that this rotating matter is 'embedded' in an inertially expanding spacetime 
coordinated frame, inertia as in Newton's first law of motion. This inertially expanding frame, however, can be associated 
with an inertial force that wants to drag (or torque) the cosmic spacetime matter out of rotation, i.e., ojfd < 0, where ujyo 
is the inertial frame dragging angular velocity in the e^-direction. It seems that the cosmic rotation is coupled with gravity 
and the inertial expansion is associated with the initial force that 'ignited' the Big Bang. 

Now, let us go further to assume that the Big Bang was perhaps due to this inertial force 'stretching' the cosmic 
matter apart as spacetime expands (like splitting the nucleus of an atom) with a cosmic cataclysmic quantum-gravitational, 
Exs = M\jC? , type explosion that caused infinitely dense matter to expand relativistically outward due to a force with a 
strength similar to that of the singularity that gives rise to the event horizon of a black hole [which draws matter inward, 
even light, due to warped spacetime according to Einstein], where M\j is the total mass of the Universe at i = 0, and E\j is 
the total energy of the Universe. It seems that this inertial force wants to 'flatten' spacetime wherein world lines will have 
straight instead of curved geodesies. This inertial force field appears to play an intrinsic part in the process from whence 
our Universe is expanded. Perhaps its origin was the initial force, as mentioned above, needed to release the infinitely large 
potential energy trapped in pre-existing Big Bang conditions in which, at least, the physical forces of nature were unified. 

The existence of this inertial field appears to be similar to the idea that Einstein had in proposing, in his Theory of 
General Relativity (see de Sitter 1916), that 'world-matter' (the matter he postulated to be the origin of the energy of inertia) 
was located at the boundary of the Universe, and controlled and provided energy for the whole Universe, from so-called 
supernatural masses. Einstein, it seems, was using his idea of conversion of matter to energy. It appears that de Sitter (1917) 
later persuaded Einstein to adopt a new hypothesis with the world-matter not at the boundary of the Universe, but distributed 



magnetism (sometimes loosely referred to as gravitomagnetism) , which is a mathematical analogy between weak gravity and Maxwell's 
equations for electromagnetism (see Mashhoon, B., 2003, [arXiv:gr-qc/0311030[ |. 
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over the whole Universe, which is finite (a sphere or an elhpse), though unhmited (i.e., no boundary). In this new hypothesis 
inertia is produced by the whole of world-matter, and gravitation is produced by local deviations from homogeneity. 

Now if we constitute Einstein's new hypothesis with properties of the proposed inertial field, then, based on this con- 
stitution, the global gravitational field of the rotating cosmic matter is a deviation from uniformity of the cosmic inertial 
expansion. If the energy of inertia is to keep the Universe expanding, then deviations from this due to gravity and rotation 
would have negligible effects locally, i.e., on small scales, but on large scales, say globally, any significant deviation may 
have an effect on the expansion rate. This has been realized in the standard cosmological Big Bang model described by the 
Friedmann- Robertson- Walker metric, and incorporated in the deceleration parameter: 

q = q(t) = qt = -^, (1) 

which depends on the Hubble constant H{t), as related to the mass (or mass-energy) density p{t) of the Universe (see 
equation I33|) . and the scale factor R and its derivatives, with A — p = 0, where, when the subscript is 0, it means t = to, the 
present epoch. 

A further example of affecting the expansion is the recently observed cosmic acceleration. Now, gravitating systems 
are associated with mass and, it seems safe to say, rotation (through conservation of angular momentum). Since the mass 
density has an effect on the spacetime expansion, thus determining the geometry in the standard cosmological model, then 
cosmic rotation (appearing to a property of gravity through conserved angular momentum) may also have an effect on the 
expansion. The apparent tendency of the inertial spacetime force (or world-matter) to stay expanding in a flat spacetime and 
the dragging of this spacetime on the rotating frame associated with gravity (the mass-energy curving spacetime) are the 
physical mechanisms proposed here to be the origin of recently observed acceleration of the comic expansion. 

Consistent with what is allowed by the Godel-Obukhov spacetime metric, the frame dragging velocity, cjfd, can be < 
or > 0, with the < solution naturally being chosen if we express lofo in the usually sense, as we shall see in the following 
section. To understand physically what is meant by ojpd < 0, we use the analogy of a rotating black hole. As mentioned 
above we assume that the Universe has two major degrees of freedom, rotation and expansion (which includes infall). Now, 
in the case of a massive rotating black hole, the gravitational force of the black hole drags inertial frames into rotation, with 
ojfo > 0. But in the case of the matter of the Universe, undergoing what one might call an 'anti-gravitational' expansion 
(i.e., a reversed gravitational collapse), the inertial force of the expansion appears to drag (or torque) the spacetime cosmic 
matter of the Universe out of rotation, with ojpd < 0. In both cases, a force is produced, which we refer to as the GM force, 
that acts on any moving matter or particle in the dragged frame. It is like say the Coriolis force acting on moving matter in a 
rotating frame, analogous to the Lorentz force acting on a charged particle in a magnetic field: hence came the word GM (see 
e.g., Thorne, Price, & Macdonald 1986). However, importantly, in the case of the Universe, with cjfd < 0, this GM force is of 
a repulsive nature, and could very well be associated with the present-day observed cosmic acceleration. Note, in the case of 
the black hole, with ojfd > 0, the GM force is also of a repulsive nature (see Williams 2004, 2002) 

The mass density p{t) regulates the warping of spacetime, while supplying the inhomogeneity to the original homo- 
geneously expanding Universe. It appears that the cosmic inertial expansion is a flat spacetime phenomenon like light (or 
electromagnetic radiation) traveling in an inertial frame, i.e., without the influence of gravity. The expansion, it seems, behaves 
like light, with general relativistic gravitation properties as well. In other words, the expansion appears to have the inertial 
properties of electricity and spacetime properties of gravitation. This would be expected, it seems, of a large scale, once unifled 
force field, of, at least, gravity and electromagnetism. 

Now, consistent with the degrees of freedom stated above, we will assume that at t ~ 10~^^ s, at least three large- 
scale forces, producing the following magnitudes of acceleration, were present to act on the mass-energy of the Universe: 
(1) the GE acceleration of gravity gcE associated with p{t), (2) the GM acceleration of gravity gcM associated with cosmic 
rotational velocity cJrot, and (3) the initial acceleration of the expansion ai associated with the inertial spacetime expanding 
coordinate frame (indicated by subscript /). Next, we consider the following scenario: If we assume that the energy, or the 
work done by the force, of the Big Bang at t = was enough to overcome the singularity binding force associated with 
qge (due to the mass-energy that initially warped spacetime closed), then we will have expansion with ai ^ (?ge, implying 
a flat or open universe. We assume that at the event of the Big Bang the Universe became gravitationally unbound, with 
matter transforming to relativistic particle expansional and rotational energy. The Universe, in general, will expand with the 
expansion velocity given to it by the force of expansion, Fj — Fi{t), which appears to be expressed by Fi ~ Euc^^H^r, 
where, again, H = H{t) is the Hubble constant, and r = r{t) is the spacetime separation between events (equation [37|. Since 
in the Big Bang 'explosion' matter was converted entirely into energy, then the expansion velocity vi, at f ~ 0, was at least 
« c, the speed of light. The acceleration of the expansion, ai ~ H^r, will decrease over time. If the scenario ended here, add 
inflation, and exclude rotation, this would be a universe explained by the standard FLRW cosmology. But with the existence 
of Qgm, of a flctitious-like force, associated with the non- inertial rotating frame, the expanding mass-energy of the Universe 
will experience an additional acceleration, perhaps like that of the recently observed cosmic acceleration. Now, this is where 
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the FLRW cosmology develops the well known problems pointed out in Section[T] i.e., when attempting to explain the physics 
of the accelerated expansion (or cosmic acceleration) we observe to exist in the present-day Universe. 

It appears that the inconsistencies in the standard cosmological model of the Friedmann-Robertson- Walker spacetime 
metric are due to our lack of considering the effects of cosmic rotation, which requires the use of a rotating and expanding 
cosmological spacetime metric, like that employed in this present paper. In the following sections the physics we need to 
further discuss the model described above, and to test its validity with observations, is devised. 



3 FORMALISM 

3.1 The Godel-Obukhov Spacetime Metric in Spherical Coordinates 

The Godel-Obukhov (Jain et al. 2007; Obukhov 2000; Obukhov, Korotky, & Hehl 1997) shear free and spatially homogeneous 
spacetime metric, defining separations of events in Cartesian comoving coordinates, is given by (Jain et al. 2007; Obukhov 
2000) 

ds^ = dt^-2V^i?(t)e'""dMy-i?^(t)(da;^ + fce2™"dy2+d2^), (2) 

where (j{t), m, and k are related geometrical parameters; R — R(t) is a time dependent scale factor, and A; > ensures absence 
of closed timelike curves (note, k is not the curvature index); with c = 1 unless noted otherwise. Clearly, ait) must be > 0, and 
for definiteness, we choose m > (Obukhov 2000). According to equation ([2]) the Universe is spatially homogeneous, rotating, 
and expanding. Note, the usual Godel (1949) metric that suffers from the presence of closed timelike curves is obtained by 
setting 

i?(t) = 1, (T(t) = 1, m = 1, fc = -^ 

in equation ([2]). The magnitude of the global cosmic rotational velocity u-^ot oriented along the z-axis is (Obukhov 2000; Jain 
et al. 2007)) 



"-' = V^^=2^\/fcT^^° (3) 



m = 2i?LJrot Y (4) 

(see equation[2}, where, recall, R — R{t) and a = o-{t). Thus, we see that, vanishing of m and/or (j{t) yields zero vorticity. 

Upon assuming a spinning fluid of proper angular momentum along the e^-axis with electromagnetic dynamical character- 
istics in a Riemann-Cartan spacetime (Obukhov 2000; Minkevich 2006), Obukhov (2000) gives an exact solution to Einstein's 
field equations: an equation of motion describing the evolution of the scale factor R. From Obukhov (2000) we can show that 



R - +3f^('^ + '^)(3. + 4.) + -^^^j(4A^-A0^ + ^^^j^c^p-p-^J, (5) 

where, the variables H , R, uj^ot, B (the cosmic magnetic field strength), p, and p, are all functions of time; Ai and A3 are 
cosmological coupling constants of the curvature, spin, and torsion tensors (Obukhov 2000) mentioned in Section[T] It appears 
that the parameters a and k, in a sense, determine the magnitude of acceleration of a fluid element due to rotation of the 
Universe; we shall see more evidence of this in Section [3. 2 1 Note, B is related to the spin density (angular momentum per unit 
volume), as we shall see in Section [5.51 Thus, one can see the repulsive nature of the above equation of motion, for the scale 
factor R, in a Riemann-Cartan spacetime, as found, it appears, independently by Obukhov (2000) and Minkevich (2005). 
Minkevich, Garkun, & Kudin (2007) have found this repulsive characteristic not only in the extreme conditions of the early 
Universe, but also at sufficiently small energy densities of later times. Minkevich et al. (2007) conclude that the effect of the 
accelerated cosmological expansion, even of today, is geometrical in nature and is connected with the geometrical structure of 
spacetime. Indeed this might be the case: it appears that these authors are finding the effect that frame dragging, producing 
the CM field, has on the geometry of spacetime. In this present paper, with the author's model proposed independently and 
unaware of concluding remark by Minkevich et al. (2007), it is shown that the recently observed cosmic acceleration may 
be the effect of the frame dragging nature of cosmic rotation, interacting with an inertially expanding spacetime geometry. 
Since the GM field is an inherent gravitational property of a rotating spacetime, it seems, we can safely say that the repulsive 
nature found by at least one of these authors (Obukhov 2000), as in the spacetime evolution of equation ((5]), is due to the 
presence of the GM field, more or less. Note, in deriving equation ([5}, for a specific epoch time t, 

H = H{t) ^H, = ^, (6) 
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the Hubble constant, was used. In addition, one cannot help but notice that the first term on the right-hand side of equation ([5]) 
is the same as that in the standard FLRW model: wherein the equation of motion of the scale factor reduces exactly to this 
term when A = p = K (curvature) = and q = 1 (see equations 1321 and I33|) . We will return to this and similar comparisons 
in Section [5.41 

For the Godel-type universe of equation ((2|, the evolution of the scale factor reveals several possible stages of the Universe 
as pointed out by Obukhov (2000), and elaborated on here, in this present paper, based on equation ([5]). The first stage is 
short and occurs in the vicinity t = 0. There is no initial cosmological singularity due to the dominating spin contribution, 
a characteristic of Einstein's gravitational theory in Riemann-Cartan spacetime (Minkevich 2006), in which R{t = 0) 7^ 
implies a regular, as opposed to a singular, spacetime metric in the transition from compression to cosmological expansion. 
The duration of this first stage is ^ 1 s, since the spin term quickly decreases with the growth of the scale factor (Obukhov 
2000), as can be seen somewhat in the third term on the right-hand side of equation But, importantly, notice that cjrot in 
the denominator will cause this accelerating term to increase over time in some degree as uj^ot — > 0. Now, at this stage (^ 1 s) 
the fluid source describing the material of the Universe is characterized by the stiff matter equation of state: p = c?p, where p 
is the pressure and p is the mass-energy density of matter and radiation. But before this stage, at very small R{t) , perhaps at 
t ^ 10~^^ s, with the equation of state possibly being that of a gravitating false vacuum (p = —(? p), however not a necessity, 
the right-hand side of equation ([5} shows that the the forth term (which appears to be related to the inertial spacetime 
expansion) together with the second and third terms (which appear to be related to cosmic rotation) will produce perhaps 
a large acceleration, one that might be associated with the inflationary stage of the standard cosmological model. During 
inflation it is commonly accepted that R{t) ~ e^*, assuming J? to be a constant during inflation (i.e., since it is approximately 
constant over the very short timescale of inflation). Next comes the stage when the scale factor increases like R{t) oc t^^^ , 
while the equation of state is of the radiation type, p « (? pjZ. This 'hot universe' expansion lasts until the Universe becomes 
matter dominated. After this the 'modern' stage starts with the effective dust equation of state p ~ 0. The scale factor still 
increases, now like -R(t) oc t^^'"'' ^ but the deceleration of the expansion takes place. The final stage depends on the value of the 
the third term referred to as the cosmological term by Obukhov (2000), containing Ai and A3, which specifically are made up 
of coupling constants relating torsion and curvature, where torsion can either accelerate the expansion or prevent cosmological 
collapse. That is, either the future evolution enters the eternal de Sitter type accelerated expansion, or expansion ends and 
a contraction phase starts. Notice the striking similarity of this cosmological term and accelerated expansion to the popular 
view of the cosmological constant A as the source of the present-day observed accelerated expansion. 

It appears that torsion of spacetime might be directly related to inertial frame dragging, and, thus the GM force field. In 
support of this, the characteristics of torsion given by Mao et al. (2007), that a rotating body also generates torsion through its 
rotational angular momentum, and the torsion in turn affects the motion of spinning objects such as gyroscopes, are exactly 
those of the GM field (see e.g., Williams 2002; Thorne et al. 1986). Therefore, it seems reasonable to associate the dominant 
repulsive nature of equation ([5]) with that of the GM field. In this paper, we derive the GM field associated with cosmic rotation 
to see what role it may have in the recently observed cosmological accelerated expansion. Further analysis of equation ([5} will 
allow us to identify, as we shall see see in Section [5.41 the suspected GM acceleration and other terms one would expect to be 
measured by a rotating and expanding comoving frame observer. 

The above stages are consistent with the model description proposed in Section (2] which includes being consistent with a 
Big Bang cosmology: first proposed by Lemaitre (1931a). This could mean that if expansion is part of the conditions occurring 
around t = lO"**^ s, then rotation, which appears to be a natural phenomenon associated with gravitation, could very well 
be a part also. So avoiding the initial singularity that exists at R(t=0)=0 for the standard FLRW cosmology, suggests that 
the Godel-Obukhov metric allows us to get somewhat closer to conditions existing at t = 0. That is, perhaps cosmic rotation 
and cosmic expansion are intrinsic parts left over from the earlier quantum-gravitational spacetime makeup of the primordial 
matter of the Universe at t ~ 0. Image that cosmic expansion and deceleration of cosmic rotation of the Universe are a 
reversed process of gravitational contraction (or collapse) and conservation of angular momentum. This helps one to conceive 
the strong possibility of how the two: rotation and expansion, cannot, it appears, be separated in the physics to describe the 
Universe, as commonly done by assuming cjiot = 0. 

The Godel-Obukhov metric of equation ([2]) converted into spherical coordinates for non-stationary and axisymmetrical 
characteristics of spacetime, which require that the metric coefficients be independent of the azimuthal (j) coordinate: g^i, = 
g^,.{t,r,e), is given by 

ds' = dt^-2y^i?(t)e'"''"°''rsinedfd<j!.-ii^(t)(dr' + r'de^-ffce'''"''='"V'sin^6'd0'), (7) 

where the usual transformation equations from Cartesian to spherical polar coordinates have been used. We identify the 
following metric coefficients for our convenience: 

gtt = 1, 

gtcf, = (7{t)R(t)e rsmO = g^t, 

grr = -R^if), 
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gee = -R\t)r^, 

94,^ = -i?'(t)fce''"''='"V'sin'6l. (8) 
Further, for our convenience, we find the corresponding inverse metric components: 

tt k 



a 



9 = 



09 

9 



1 



RHty 
1 



i?2(t)r2' 

1 



i?2(i)e2™-=i"er-2 sin^ 6l[fc + a{t)] ' 
where, in general, from inverse matrix operations, 

9il><l> 



(9) 



9 

tc, 

9 

r'l 

9 



9ttg^4- - 9% ' 



gt4> _ ipt 



9tt94>4> - 9^4, 



1 

Qrr 



9 



1 

909 ' 

</)<p _ 9tt 



99 

9 



-■ (10) 

9ttg4,4, - 9t4, 

The frame dragging velocity for the Godel-Obukhov metric of equation 0, with diagonal component signature (+, -, -, 
-, -), is given by 

9't't 

CJPD — 

9<l'<t> 



i?(t)A;e™''"ner sine' ^^^"^ 

where we are being consistent with the geometrical expression given by Bardeen, Press, & Teukolsky (1972) for the frame 
dragging velocity. In equation we see that frame dragging velocity in the direction can be either positive or negative 
because of the square root, but cjpd < occurs naturally it seems, producing, as we shall see in Section 13.21 a repulsive 
gravitational force. 



3.2 The Cosmic 'Gravitomagnetic' (GM) Force 



We now derive an expression for the GM force, J^'gm , exerted on a test particle (or an object such as a galaxy) of cosmic space 
momentum P. Using the analogy of expressions found in Williams (2002, 2005), for a rotating general relativistic system such 
as our Universe, the GM force measured by a comoving observer can be given by 





dp 




( 


dr 


) 






GM/ 



— H 



dP 

d7 



HP 



(12) 



(dr is the proper time interval), with 

Hrj = e""(^GM),|. 

(the vertical line indicates the covariant derivative in 3-dimensional absolute space), where, like in the case of a rotating 
black hole (Thorne et al. 1986), 

(/3gm)'' = (/3gm)'=0, (/3<,,J^ = -^fd; (13) 

u)-po is given by equation pip . The field H is called the GM tensor field, and fi^^ is sometimes called the GM potential. We 
perform the metric component operations in eq. p2p to give the following expression for the GM force exerted, Pqm: 



dP 

dr 



= Fc 



— [(-fGM)r, (-fGM)e, (-fGM)^] 
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(14) 



Here we are only interested in the radial component of eq. (|14p . where we are assuming that the other components are not 
important in explaining cosmological acceleration along the line-of-sight of the observer. Therefore, we need to determine the 
force in the radial direction: 



(FGM)r = HreP^ + Hr4,P'^ 

= Hreg^'Pe + Hr^g'^'^P^ 



(15) 



where we have used P^ = g^" P,,. (Note, equation [14] is a general expression, existing for any gravitating and rotating system.) 

In equation (|12p we can identify the so-called blueshift factor e~" (= \fg^) for a metric with a signature of diagonal 
components of the type as in equation ([7]). Thus, from equation Q, 



k 



k + a{t) ' 

Next we determine the relevant GM tensor components: Hre and H^^, to be substituted into equation (|15|l . These 
components are given in equation p2p : 

(17) 

where we have used equation (|16p . In general the covariant derivative in 



(repeated indices of i, j, k sum over r, 
3-dimensional absolute space is given by 

SO that 

Hre = y5"(/3GM)e|r 



= -Vff"rer(/?GM)0, 

since, as given by equation (|13|l . (/3GM)r 

(/3gm),;> — gii>4, (Pgm)'^ 
= g^0(— ojfd) 
= =Fy^P(t)e'"''""Vsin( 



(18) 



(19) 



(;5GM)e = and (^gm)^ / 0, where 



(20) 

upon substitutions from equations ((8| and (|11[) . Note, (/3gm)(^> is measured in units of length and (/3gm)'* in units of per length. 
Note also that for lufu > 0, (^Sgm)"^ < (= GM potential), and (/^gm)^ > 0; for ujfu < 0, (/Sgm)** > 0, and (Pgm)^ < 0. Now 
similarly. 



■ (/3GM)</,|r 

^ [(/3GM)0,r - r^^(/?GM)0] 



9(/3gm),^> 
9r 



r^,,(/?GM)0 



(21) 



It can be shown that Hre ~ 0, from equations ((HJ, ((9]), ([19]), and ([24]); therefore, the GM force in the radial direction 
relative to an arbitrary comoving observer, given by equation (|15[) . reduces to 



(PGM)r = -ffr^S -P.* 



(22) 



We want to simplify and analyze the above vector component to see under what conditions it may contribute to a repulsive 
accelerating force, i.e., we want the GM radial force component to be repulsive (> 0). We first evaluate the partial derivative 
of equation (|2H) . The negative value of equation (|20|) is chosen, consistent with the note stated above following equation (|20[) 
concerning ojfd < 0, and, importantly, consistent with the description of the model presented in Section [2] concerning the 
frame dragging, ljfd < 0, which tends to drag the Universe out of rotation. Thus, from equation (|20p . 



^(^gm)^ ri7\ Du\ ■ a ^ / 

= -y'a{t)R(t)sme—{re 



mr sin ( 



/a{t)R(t) sine e 



mr sm ( 



(mr sin 6 + 1). 



We next evaluate F^^ of equation (|21|l . In general the afhne connection is given by 



1 X. 
o5 



dgni, ^ dg^ 



(23) 



(24) 



Could Dark Energy he a Manifestation of Gravity? 9 



So, with gtr = ff0r = 3 



1 



Ipt 



dgt. 



2^^ y dr 
mr sin 6^ + 1 



4,4, 



2r 



+ 2^ 
'cr(t) + 2k 



dr 



k + a{t) 



upon substitution of nonzero metric components from equations 
equations ^3^, ((25l), and ([20} into equation ((2l| yielding 



(25) 



and Q into equation (|24p . Now we substitute from 



ka{t) 



k + a{t) 



1/2 



R{t) sin 9 e 



mr sm t 



(mr sin 61 + 1) 



a{t) + 2k 
2[k + a{t)] 



- 1 



(26) 



Finally, upon substitution of equations (|26|) and ((9]) into equation (1221) yields the following for the cosmic GM radial force 
component along the line-of-sight: 

1/2 



(Fa 



ka{t) 



[k + a(t)]3 



mr sin 9 + 1 



R(t)r'^ sin Oe 



mr sm 6 



a{t) + 2k 
2[k + a{t)] 



(27) 



in geometrical units (G = c = 1), where cr{t), R{t), and k are dimensionless; and m has unit of per length, with the scale factor 
being normalized at the present epoch: R(t = to) = 1 (see equation I43|) . Since the angular momentum of a cosmic object of 



mass M and cosmic circular velocity v^, is P4, 

1/2 

I tcm T 1 I 

(-FGM)r = ± 



L — Mujrotr , then, expressed in proper units for calculation purposes. 



fecr(t) 



[k + a{t)]3 



{m/c)r 



R{t) sin 6»e('"/<=)'-=in9 



1 - 



a{t) + 2k 



Mcuiotc, 



(28) 



2[fc + a(t)] 

where, again, ojrot is the global cosmic angular (or rotational) velocity, and m is given by equation The positive of the 
above expression is for cjfd < 0, and the negative is for lufu > 0. According to equation (|lip . it appears that only cjfd < 
is physically sound, i.e., (Fgm)i- > 0. 



3.3 The Cosmic 'Gravitoelectric' (GE) Force 

We now calculate the cosmic gravitational force throughout an assumed spherical axisymmetric universe (of infinite extent 
to an arbitrary comoving observer). We will assume that the equation of state and scale factor relations, derived from the 
Einstein field equations, by way of a FLRW cosmological model, are still at least approximately valid in the Godel-Obukhov 
cosmology. Support of this assumption is that the rotating and expanding Godel-Obukhov metric (equation [7|, in the limit 
of large times and nearby distances, reduces to the open metric of Friedmann (Carneiro 2002). Moreover, we will also assume 
that the derivation of the GE force for the spherical axisymmetric case is not much different from that of the spherical 
symmetric (FLRW) case. The result will allow us to test the validity of this assumption once the exact GE term is identified 
in equation Q. 

Note, the Godel-Obukhov metric is shear free, spatially homogeneous, and isotropic in the cosmic microwave background 
(CMB) radiation (like in the standard FLRW cosmology) for any moment of cosmological time t (Obukhov 2000). Further, 
and in summary, the cosmological model of equation 0, with rotation and expansion, does not suffer from the three major 
problems associated in the past with cosmic rotation. This cosmological model is causal, isotropic in the CMB radiation, 
and parallax free; and thus, the limits on the cosmic rotation, obtained earlier from the study of CMB radiation and of the 
parallaxes in a rotating world, are not true for the class of cosmologies in which the Godel-Obukhov metric of equation ((Tj is 
a member (Obukhov 2000). 

The gravitational potential inside the Universe is assumed to be given by the post-Newtonian approximation 

$(r) = -G / dV ^ ' (29) 
J |r — r 

where the integral is over the distance from the arbitrary observer located at r' — robs = 0, along the line-of-sight, out to an 
arbitrary distance r; and where 

T°° ^J2m„5^r-r') (30) 

n 

for a gravitational bound system of masses m„. The component T"" is the rest-mass density, or commonly referred to as the 
mass density, of the energy-momentum tensor, y^"^ which serves as the source of the gravitational field. For nonrelativistic 
matter equation (|30|l can be set equal to the mass density p{r'). Again, we are assuming spatial homogeneity, and isotropy 
in CMB radiation only, which we will refer to as the 'quasi-Cosmological Principle.' In the FLRW model, the Cosmological 
Principle of spatial homogeneity and spatial isotropy is assumed, which is consistent with CMB temperature measurements 
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[save for the puzzling anomalies found in the WMAP temperature maps that are not expected from gaussian fluctuations 
(Chiang et al. 2003; de Oliveira-Costa 2004)], and consistent with large-scale structure observations [save for the large-scale 
asymmetries that are equally unexpected in an isotropic, homogeneous space (Eriken 2004; Land & Magueijo 2005)]. These 
measurements and observations confirm however, to a strong degree, the Cosmological Principle. Perhaps the small anomalies 
and asymmetries are effects predicted by cosmic rotation that do not conform to the standard FLRW cosmological model. 

For a universe where the centre is everywhere, according to the Big Bang, it seems reasonable to assume that at any 
given time the Universe can be represented as a rotating sphere of average constant density pt, such that for spacetime 
expanding from a Big Bang origin, the distance the Universe has expanded from its initial 'point' (i.e., state) is equal to 
the coordinate separation between galaxies (or protogalaxies). The above reasoning allows use of equation 1)29^ to derive the 
GE (i.e., the familiar gravitational) acceleration experienced by all galaxies throughout a spacetime expanding universe. The 
requirement must be that for a freely falling, locally flat spacetime observer this gravitational acceleration is approximately 
zero, according to the Equivalence Principle. So, evaluation of equation (|29|) . assumed to be that measured by a comoving 
observer at a coordinate separation distance r, with a constant mass density throughout the Universe, p = p(t) = pt, for any 
given expanding epoch, yields the following expression for the universal gravitational force exerted on a mass M at a distance 
r from an arbitrary observer at robs as defined above: 

(Fge). ~ - I^GprM, (31) 

where by deductive reasoning r is a measure of spacetime separation, identified as the same r as in the spacetime metric of 
equation ((T]). Dividing equation (|31|) by M gives the universal gravitational acceleration for a given epoch at a distance r 
from an arbitrary comoving observer. Equation H31[) satisfies the above requirement that the gravitational acceleration goes 
to zero as r — >■ 0, as measured by the comoving local inertial spacetime observer. The validation of the above reasoning used 
in deriving equation (|31|) will be given in the following section. Importantly, we shall see that pt is just the mass density in 
the equation of motion of the scale factor in the standard FLRW cosmology when A = = p = (see equation I32[). Notice 
that the GE force given by equation (|3ip is negative and opposite the sign of the radial component of the GM force given 
by equation ((28]). This repulsive nature of the GM force is consistent with it acting to accelerate the cosmic expansion of the 
Universe. To test this claim of consistency, in Sections [4] and O we will compare the magnitudes of the accelerations produced 
by the GM and GE forces at redshift z ~ 0.5, to see which is dominant. 

To summarize, we are assuming that the Universe can be represented by a spherical axisymmetric cosmology (though 
infinite in extent to a comoving observer). The force i^cE of equation (|3ip expresses the gravitational force, due to the average 
mass density pt (being a constant throughout the Universe for any given epoch according to the quasi-Cosmological Principle), 
acting on say a galaxy of mass M at a distant r, as measured by an arbitrary observer, where, for this observer, r — > 0, which 
means that Fqe — > 0, as it should locally, in accordance with the Equivalence Principle, assuming any arbitrary comoving 
observer in the Universe is in free fall. Notice, however, the same is not true for (Fgm)i- of equations (|27p and (|28p . i.e., f gm 
does not go to zero at the observer, where r — > 0, because Fqm exerts a force on locally moving inertial frames; then only 
if oirot — > will Fqm — > 0. In other words, the GM force in general acts on the momentum of a test body or particle 
[with effects being readily seen in a strong gravitational potential such as that of a black hole (e.g., see Williams 2004)]. The 
findings in this present paper suggest that the GM acceleration has a significant dynamical effects on spacetime comoving 
objects in a rotating frame, and, thus, on the expansion rate of the Universe. 



3.4 The Density of the Universe 



We now derive an expression for the mass density pit) of the Universe, which includes any contribution from radiation. 
We assume that the standard FLRW cosmological model is approximately correct. The Friedmann-Lemaitre's solutions to 
Einstein's gravitational field equations yield the following acceleration equation: 



R 
R 



A 
3 



4nG 



where the Robertson- Walker metric was used and K {= 
critical mass density is given by 



Pc{t) 



47rG ' 



(32) 



curvature) = 0. Then the general expression for the time-dependent 



(33) 



with A = p = 0, implying specifically a Friedmann cosmology, where we have used equations fTJ and ((6|; pc is the density 
needed to make the Universe flat. Note, the way in which pc of equation (|33p was derived, from the standard FLRW cosmology 
(equation I32[). does not rule out contribution from dark energy, but only sets A = 0, for an isotropic universe: no pressure 
gradients (Vp — 0) and nonrelativistic matter or matter dominated (p — 0) universe. This does, however, suggests that dark 
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energy might be directly related to gravity and that the standard FLRW cosmological model cannot fully account for its 
presence. Thus, the critical mass density pc in terms of the measured present epoch cosmological parameters, is given by 

Pit = to)c = Pc(to) = ^g^, (34) 
where we find that 

Pc(to) « 9.5 X 10"^° gcm"^ (35) 

for the currently suggested values of Ho ~ 71 kms~^ Mpc~^ and go — 5. This value of deceleration parameter qo indicates a 
flat universe, which implies that Q, = p{to) / Pc{to) — 1, consistent with observational data finding that f2 ~ 1 (Jarosik et al. 
2011; Tegmark et al. 2004), if we assume that p = pmat + pa; then division by pc yields Q = flmat + ^lA — 1, where observations 
suggest that fimat — 0.27 and f^A — 0.73, requiring A 7^ and p < in the standard FLRW cosmology, specifically defining 
the ACDM cosmological model. This model, however, fails to tell the true nature of so-called dark energy, leaving the subject 
open to speculation. Nevertheless, equation (|33[) can be identified as the source of the universal gravitational field of attraction 
and will be used in the GE force given by equation (|3ip . which gives the GE force for different epochs, where evaluated using 
equation (|34|l gives the present strength. 

Now we return to give validity to the expression for the universal gravitational force of attraction (equation I3ip . and, 
thus, to the reasoning that led to its derivation. Upon substitution of the critical density of the Universe (equation I33[l 
into equation (|3ip and dividing by the mass M, we obtain the GE acceleration [i.e., the gravitational force per unit mass 
(gGE)r — r]. We can express this as a deceleration of the cosmic scale factor R: 

r ~ —qH^r, 
RX ~ -qH'Rx, 

f ^ -qH\ (36) 

which is, as would be expected, the same as that of the standard FLRW cosmology, when A — p = and equation (|33p is 
substituted into equation (|32p . where 

r{t) = R{t)x, (37) 

relating the physical distance r to the comoving coordinate distance Xi s^^id its derivatives, have been used. The vector x 
comoves with the cosmic expansion. One can think of equation (|37p as a coordinate grid which expands with time. Galaxies 
remain at fixed locations in the x coordinate system. The scale factor R{t) then tells how physical separations are growing 
with time, since the coordinate distances x s-re by definition fixed. Further, solving for q, we can identify equation (|36|) as 
that of equation ((T]), with H given by equation i.e., we identity the deceleration parameter as defined in the standard 
FLRW model. Again, this is what one would expect for the behavior of the GE force of equation ((31]), as it relates to the 
standard model, and, thus, this can serve to validate the reasoning behind assumptions made in its derivation. Importantly, 
note, equation (|36p is exactly equal to the first term on the right-hand side of equation ([5]) with g = 1. Therefore, this term 
can be identified as the GE deceleration of the scale factor in the Godel-Obukhov spacetime (we will return to this discussion 
in Section O. So, in summary, the validity of the derivation leading to equation pip , which can be used to express the GE 
acceleration approximately in both the Godel-Obukhov and FLRW cosmologies, has been established. The assumption that 
the derivation of the GE force for the spherical axisymmetric case is not much different from that of the spherical symmetric 
(FLRW) case has been validated, at least qualitatively; and it appears from equations ((5)1 and (I36p that the strengths will 
differ quantitatively by a factor of q. 



3.5 Cosmological Parameters 

For a qualitative and somewhat quantitative analysis of the model described in this paper we choose the following parameters 
of equations ([Sjl, (|4)|, and ((T}: a, m, k, and tJrot, based on observations and, of course, on theoretical insight. A possible way 
to express evolution of the force Fgm of equation (|28p over time is to let 



cr = a(t) EEe"i*/*°, (38) 
and let k be defined as a function of a{t) by Obukhov's (2000) model relation 

k = cia{t), (39) 

where, when estimated from the Godel-Obukhov metric and general relativity, C2 ~ 71, using go = 0.01, (ajrot)o = O.li^o, 
Hi) ^ 70 kms~^ Mpc~^ (again, the subscripts indicate the present epoch), with tJrot = ajrot(t) = {y>ro\)t- A reasonable choice 
for the constant ci appears to be ci ~ —115. The value of ci is related to the magnitude of the force f gm of equation (|28p : 
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for example, upon changing from ci — —105 to ci = —115, in these model calculations, the magnitude of the force, for a 
typical case, increases by about two orders of magnitude. Note, with such choices, at f = 0, cr(t) ~ 1 and fc ~ 71, consistent 
with the fe ^ requirement for causality. We will use the value C2 ~ 71, since making it relatively larger of smaller appears 
to have little effect on the model outcome. 

We will assume an analytical expression f equation 154 |l consistent with the more recent estimate, given below, for the ratio 
of the magnitude of cosmic rotation to the Hubble constant, where observations of anisotropy in electromagnetic propagation 
from distant radio sources, expected typically of cosmic rotation, are used to determine the estimate (Obukhov 2000; Obukhov 
et al. 1997; Nodland & Ralston 1997a): 

%i^=6.5±0.5, (40) 
with galactic coordinate direction I = 50° ± 20°, h — —30° ± 25°. This value is larger than the estimate (Obukhov 1992): 

%i^ = 1.8±0.8, (41) 
Ho 

with direction / = 295° ±25°, b = 24° ± 20°, obtained from Birch's (1982, 1983) data. Moreover, recent analysis of the 
large-scale distribution of galaxies (Broadhurst, Ellis, Koo, & Szalay 1990) has revealed an apparently periodic structure 
of the number of sources as a function of red shift. From this we get yet another estimate of the rotational velocity which 
appears necessary to produced this observed periodicity effect. This estimate gives the largest ratio of the three (Obukhov 
2000; Korotky & Obukhov 1994): 

^ 74 (42) 

(compare equations l40l and l4H) . It is clear from above that further careful observations and statistical analyzes will be extremely 
important in overcoming the inconsistencies, in establishing the true value of the cosmic rotation (or vorticity), which may 
result from too few empirical data. 

Next, we use our defined scale factor R{t) to relate the Hubble constant with time. In general, with assuming the usually 
power-law solution for the scale factor as a function of time {R oc t") according to the FLRW cosmological model, based on 
the the assumptions of Section 13.31 




Rit) = j , (43) 
normalized at the present epoch t = to. Using equations ([l]) and ([6|, we get the general expressions 

Ht = n't-\ (44) 



and 



(45) 



where, in equation (I44|) . n = n' , which gives for n' = n = | an age of the Universe (9.22 x 10^ yr) too low to be consistent 
with recent observational estimates of Ho — 71 kms~^ Mpc^^, with <Z = | according to equation (|45p . On the other hand, the 
expression Ht = t^^ , n' = 1, gives an age (13.8 x 10^ yr), which is consistent with recent observational estimate with cosmic 
acceleration (Lineweaver 1999), and without acceleration, in the absence of deceleration (Freedman et al. 2001). So it seems 
reasonable to assume the following limits for the present age to: 

^<io<;^, (46) 



3 



1. Note, with n = 1, then according to equation (|45|l q = 0, implying an open universe in the standard FLRW 
cosmology (compare equation [T|. Specifically, for concreteness, it appears appropriate to choose to = H^^ (n' — 1) for the 
present epoch, but with n = | in equation (|43p . Note, with observations suggesting that the age of the Universe is closer to 
the Hubble time (Hg^), instead of that given by the standard FLRW cosmological model (|_ff^^), implies that the Universe 
has at least not decelerated continuously. The discrepancies leading to the limits above can possibly be attributed to the 
evolution of R{t), i.e., how it might change as the Universe undergoes phase changes, thus reflecting how the value of n might 
change, where n = |, recall, is also the starting scale factor at t ~ of Lemaitre's (1931b) expanding cosmological model. 
Moreover, for completion, during inflation (indicated by the subscript 'infl') it appears that 

= e^ / ~ e^'"''"^* (47) 

[it(t)Jinitial 

(i.e., e/'^''* « e"!""'), where [R{t)]inMa.i and [-R(t)]finai are the initial and final scale factors before and after inflation; 
Hina — 1/finfl is the Hubble parameter at the onset of inflation, which remains approximately constant during inflation; and 



Could Dark Energy be a Manifestation of Gravity? 13 




1e-40 1e-30 1e-20 1e-10 
t (seconds) 



1e+10 



1e-08 1e-07 1e-06 1e-05 1e-04 0.001 0.01 0.1 

t/to 



Figure 1. The scale factor R{t) = (t/to)" versus time up to t = to = 13.8 X 10® yr, with n = 2/3 or n = 1/2. (a) Schematic plot of R(t) 
vs. t in seconds, from the Planck time (~ lO""*^ s) to to- Inflation is indicated by the steep rise in the curve at IQ-^*^ s ^ t ^ 10-3* s, 
where R{t) increases by a factor ~ e^® from its initial value given by R{t) = {t/to)", with n = 2/3 before inflation (see text), (b) 
R{t) = {t/to)" vs. t/to for 138 yr ^t ^ to- The step-like feature indicates where n = 1/2 (just after inflation) changes to n = 2/3, at the 
radiation-matter equilibrium time icq ~ 1.7 X 10^^ s (see text); this can also be seen in (a). 



At = ffinai — ^initial, wlth tinitial and tflnai indicating the beginning and ending times of inflation, respectively. For example, 
assuming that inflation occurs between 10"^'' s < f < IQ-^" s, with to = 13.8 x 10® yr, we find that the scale factor by which 
the Universe increased during inflation is [i?(t)]flnai — e®®[_R(t)]initiai — 10*3[7?(t)]initiai. Now, whether or not inflation actually 
occurred, or its origin, we do not know for certain, but we do know that the Universe, early on in its history, appears to have 
increased or inflated by a factor of ~ 10''^ from a small causally connected comoving region of spacetime ^ 3 x 10~^^ cm 
at t^ IQ-^^ s to a size ~ 3 x lO^'^ cm at t ~ IQ-^* s, according to equation p7[) and the above relation between the initial 
and final scale factors. In addition, importantly, cosmic vorticity enhances the inflation, i.e., when the vorticity is large, the 
inflation rate is much bigger than in the vorticity- free case (Obukhov et al. 2002). 



4 NUMERICAL MODEL RESULTS 

For comparison and completion, plotted in Fig. 1 is the cosmic scale factor. Figure 1(a) displays a schematic plot of the scale 
factor given by equation (|43|) . with n = | or n = |, over a period from when the age of the Universe was ~ 10"*^ s to the 
present estimated age of to — 13.8 x 10® yr. The lower time limit corresponds to a Hubble radius ~ 3 x IQ-"^"^ cm, of the 
order of the Planck scale, expanding out to a radius ~ 1.3 x 10^* cm. The Hubble radius, mentioned above, is related to 
the causally-connected observable Universe and is given by ru = c/Ht. Immediately following the Planck era we believe that 
the Universe was at least a thermal causally connected spacetime gaseous plasma. We assume that n = | in equation (|43|l 
at t = 10"*'' s, with this value lasting up to the beginning of inflation, as indicated in Fig. 1(a). Here also we are assuming 
that at the onset of inflation the Universe was at least radiation dominated (setting n = i), where at this epoch it appears 
a 'false vacuum' or the release of a type of quantized-gravity binding-like energy, resulting from symmetry braking of at 
least three of the fundamental forces (strong, gravitational, electromagnetic), drove inflation. The details as to what initiated 
inflation are yet to be understood; at present we can only speculate. After inflation, during the radiation dominated era, 
from when the age of the Universe was t ~ lO-^"* s to t = icq ~ 1.7 X 10^^ s ~ 53, 968 yr, we continue to set n — ^ 
in equation (|43|) . where tcq is the time of matter and radiation equality (Liddle 2003). Beyond tcq we set n = |, indicating 
matter dominance, producing the step-like feature clearly seen in Fig. 1(b) at t = toq ~ 4 x 10"^ to- Before this time relativistic 
particles dominated. As the Universe continues in a matter dominated phase after recombination, at t ~ 350, 000 yr after 
the Big Bang, in equation (|43p we still have n = | up to the present epoch. Note, this expression for R(t) indicates a flat, 
decelerating universe (i.e., g = | > 0, using equation I45p . as would be expected in a FLRW expanding cosmology. Yet, this is 
only somewhat consistent with observations, because recent observations appear to indicate an open, accelerating universe, at 
least for the present epoch, with g < according to the standard FLRW cosmological model. This would cause R{t) to have 
a somewhat steeper incline (or slope) near the present epoch than that displayed. For example, for q — —0.2, equation (1451) 
gives n — 1.25. Now, Fig. 1(b) displays R{t) of equation (|43p over the time (138 yr ^ t ^ 13.8 x 10^ yr) that the gravitational 
accelerations, the GM and the GE, are calculated from the forces of equations (|28p and pip , respectively, as we shall see 
below. The lower time limit is set here by the computational capacity of the computer in the units used in calculating the 
GM acceleration from equation (|28p . This limit will, however, be overcome in Section (5) Notice, the step-like feature is an 
'artifact' resulting from normalizing R{t) at the present epoch (compare equation I43p. 

Displayed in Fig. 2 are the evolutions of the magnitudes of the cosmic gravitational accelerations (force per unit mass) 
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Figure 2. The magnitudes of the radial accelerations {gGM)r (solid curve) and {gGE)r (short-dashed curve) produced by the gravito- 
magnetic (GM) and gravitoelectric (GE) forces, respectively, in cgs units, versus t/tQ, representing the evolution from t = 138 yr after 
the Big Bang up to the present time t = to: (a) Evolution of accelerations at a distance with z = 0.65 (see text). Note, the (3GM)r 
component extends down to ~ 2 X 10~^^ cms~^, then begins to turn up. (b) Evolution of accelerations at a distances with z = 0.55. 
(c) Evolution of acceleration at a distances with z = 0.5. (Note, the x-axis extends to 1.2 to so that we can clearly see the values of the 
functions at t = to, i.e., at the present epoch.) 




and 



GE)r, from equations 

1/2 



and H31|) . respectively, versus t/to, where 



(ffGM)r = 

and 

(gGE)r ^ 



ka{t) 



[k + a{t)]3 



{m/c)r sin 6 + 1 



R{t) sin 6'e(^ 



%l c)r 



a{t) + 2k 

2[k + G{t)] 



(48) 



(49) 



with 6* = ^ in equation (|48|) (see discussion in Section [5751 



The above radial gravitational accelerations, ((?GM)r and {gGE)r, 
are measured at a coordinate separation distance r, corresponding to a particular redshift z, by an arbitrary comoving observer, 
as this distance expands over time, while the gravitational accelerations at that distance evolve over time, from t ~ 138 yr 
after the Big Bang to the present estimated age of the Universe: to — 13.8 x 10*' yr. 

In these calculations we step through the independent variable t by assuming the following: 



(50) 



where fi is the fraction of the total time that we step through, normalized to equal one at the present epoch, i.e., t = to; and 
subscript i indicates a step size. The Hubble constant then evolves as 

Ho 



H{t)^ 



(51) 



The evolving mass density p — p{t) of equation (|49|) is thus given by equation (|33[) for a flat universe according to the standard 
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Figure 2 - continued (d) Evolution of accelerations at a distance with z = 0.45. (e) Evolution of accelerations at a distance with z = 0.35. 
(f) Evolution of the accelerations at z = 0, where (gGE)i- — > 0. Note, {gGM)r reaches a maximum finite magnitude at 2 = r = 0, as 
measured at the arbitrary comoving observer (see text). 



FLRW cosmology. Similarly, the evolving distance r — r{t) is assumed to be given by 

for recession velocities ^ c, which is just the nonrelativistic Hubble law. Note, for z ^ 1, the relativistically corrected Hubble 
law (equation [56]) must be used for accuracy; this will be discussed further in Section [5.21 

The above evolving distance r{t) is for a specific z, measured by a present-day observer, indicating how a specific comoving 
coordinate point in spacetime has evolved. Substitution of the evolving variables: r{t), a{t) (equation I38p . R{t) (equation I43p . 
and p{t) (equation I33|) . into {gGM)r and {gGE)r, above (equations 1481 and I49p . allows us to see how these cosmic gravitational 
accelerations have evolved at that specific comoving coordinate point, indicated by z, as measured by an arbitrary present 
epoch observer. Note, it shall be interesting to see what happens to ((?GM)r as t approaches zero and what role it may play in 
the cosmic infiationary era. It is suspected that it may play a significant role. We shall see if our suspicions are true in Section[5] 
where we will attempt to go back in time as far as theoretically possible using a valid approximation to the GM acceleration 
of equation (|48p . However, for the present, we find that for z — > in equation (|48p . {gGM)r reaches a finite maximum value of 
(<7GM)r ~ 3.9 X 10^^ cms~^ at an arbitrary comoving observer as shown in Fig. 2(f). Now, by equation (|48p and conservation 
of angular momentum, as t decreases, in Fig. 2, for t < 0.01 (or < 1.38 x 10* yr after the Big Bang), {gGM)r increases with 
increasing ojrot as would be expected. Yet, on the other hand, as t increases, for t > 0.01 to, (gGM)r first decreases as would 
be expected, but then as ojrot gets smaller and smaller, {gGM)r once again increases, at least for the z values shown (compare 
Fig. 2; see also Fig. 3). [The behavior of ((?GM)r for larger values of z will be discussed in Section [5.2j . The above behavior 
of {gGM)r appears to be consistent with the third term on the right-hand side of equation ([5]); and, importantly, in Fig. 2(c), 
the magnitude of ((?GM)r overtakes that of ((?GE)r, indicating a net positive acceleration or repulsive force per unit mass. We 
shall return to this discussion in the following section. 
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Figure 3 displays how the cosmic rotational velocity decreases over time: Specifically plotted, as we shall see below, is 
a derived analytical expression consistent with equation (|40p . We shall see that these model calculations suggest that the 
rotational velocity LJrot has a value ~ 6.3ff, which is consistent with present-day observations. In what follows, an analytical 
expression is derived for the cosmic rotational velocity cjrot of equation (|48p . whose numerical value is consistence with 
observations: As usual the magnitude of the angular velocity (or the angular frequency) is given by 

~ dt 

f , (=3, 



assuming simple harmonic-like motion. From equation (|53p . it seems reasonable to relate the Hubble constant, the rate of 
cosmic expansion, to LOrot, the rate of cosmic rotation, by 

~ 6.3H, (54) 

where we have use t ~ i. Importantly, this expression is consistence with observations; compare equation (|40p . Equation (|54p 
is plotted in Fig. 3. 



5 DISCUSSION 

In the following sections, we analyze further and discuss the results above. We will look at the long term behavior of the GM 
and GE accelerations over time. Because of the large interval of time covered, we do this in two separate epochs, to compensate 
for the limited capacity of the computer. However, since these calculations are computations of analytical equations, none of 
the physics is lost, because the values and times converge at the 'interface,' indicative of a specific cosmological time. 



5.1 The Gravitomagnetic (GM) and the Gravitoelectric (GE) Accelerations: From 138 yr to Present 

We first analyze the GM and the GE accelerations over the time for which we have the exact analytical expression for 
the GM acceleration (equation I48|) . with the GE acceleration given by equation (|49p . As expressed in Fig. 2, this time 
is 138 yr ^ t ^ 13.8 x 10^ yr after the Big Bang. Note, equation (|48p is referred to as the exact in comparison to the 
approximation to this equation we will use in the following section to find its value in the early universe (t < 138 yr). 
As mentioned in Section [l] the evolution of the GM and GE accelerations over time at a distance r{t) as measured by a 
present epoch arbitrary comoving observer for a specific z (see equation I52|) as related to the cosmological distance given 
by equation 1)37^ is plotted in Fig. 2. These model calculations show that the GM acceleration (equation [48]) goes to zero at 
2^0.7 as measured by an arbitrary comoving present observer. This implies that the Universe was in a decelerating phase 
for z at least greater than ~ 0.7, i.e., at an earlier cosmic time; this also is consistent with observations (Riess et al. 2001). 
Figures 2(a) and 2(b) seem to show that the Universe starts to decelerate at a slower and then an even slower rate at z — 0.65 
and 0.55, respectively. This is because according to Fig. 2(a), ((?GM)r begins to increase at ~ 8 x 10^ yr after the Big Bang, 
from a minimum value ~ 2 x 10~^^ cms~^, at the spacetime coordinate point, r, associated with z = 0.65 (as measured by 
a present epoch observer); and according to Fig. 2(b), ((/GM)r begins to increase at ~ 6 x 10® yr after the Big Bang, from 
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Figure 4. Magnitudes of gravitational accelerations, {gGM)r and (gGE)r, over time: (a) From the Planck time [t ~ 5 X 10 s up to 
t ~ 4 X 10^ s (= 138 yr)]. Solid curve is for (gGM)r, n = 2/3 implying q = 1/2 (see text); dotted curve is for (gGM)r', n = 1/2 implying 
q = I; dashed curve is for {gGE)r, n = 2/3; dash-dotted curve is for (gGE)r, n = 1/2. Note, (gGM)r, n = 2/3 and {gGM)r, n = 1/2 
converge at the present epoch, to ~ 4.3 X lO^'^ s ~ 13.8 X lO'* yr, since ((/GM)r depends on the scale factor, which is normalized at t = to. 
There is a factor of 1/2 difference between (gGE)r, " = 2/3 (implying q = 1/2) and {gGE)r, n = 1/2 (implying q = 1), where n = 2/3 
is for matter dominance and n = 1/2 is for radiation dominance, (b) From t = 10~* to (= 138 yr) to to = 13.8 X 10^ yr. Solid curve is 
for (5GM)r, short-dashed curve for (gGE)r, evaluated at the Hubble radius, being consistent with (a). The step-like feature on the curve 



for (gGM)r indicates the change from radiation dominance to matter dominance at t 
n = 2/3. Note, (gGM)r — > at t ~ 0.6 to (not shown; see text). 



toq ~ 54, 000 yr (see text). 



1/2 goes to 



a minimum value ~ 9 x 10~^^ cms~'^, at the spacetime coordinate point associated with z = 0.55. Importantly, we see that 
Fig. 2(c) is consistent with recent observations that suggest that the Universe entered into an accelerating phase at z ~ 0.5, 
with ((7GM)r > |(ffGE)i |. Figure 2(d) shows how as z gets smaller, which means that the distance from the arbitrary comoving 
observer is getting smaller, the GM acceleration gets larger and larger; and we find that, as z — > 0, the GM acceleration will 
continue to get larger until a maximum value of {gGM)r ~ 4 x lO'^* cms~^ is reached at the comoving observer, as mentioned 
in Section |4l and as can be seen in Fig. 2(f). 

Now looking at Fig. 2(c) in more detail and using what is known from particle physics, up to t ~ 2.7 x 10"^ to, where, 
again, to = 13.8 x 10^ yr, the Universe was filled with ionized gas. Then around this time, t, recombination occurred with 
the decoupling of radiation and matter, and the hydrogen atom formed causing the Universe to become transparent while 
entering into a 'dark-age' phase. The fundamental force of gravitation became dominant on small scales, eventually leading to 
large-scale structures, as the Universe continued to cool. The so-called dark ages persisted up to around t ~ 4 x 10~^ to, until 
quasars and/or protogalaxies formed, and then the Universe under went a 'cosmic renaissance' phase, ending the dark ages, 
and entering into a reionization (of hydrogen) phase. Figure 2(c) shows that the magnitude of the GM acceleration becomes 
less than that of the GE acceleration during the so-called dark ages at t ~ 1 x 10~* to, indicating that the Universe goes into a 
decelerating phase at the evolving cosmological distance r{t), associated with z = 0.5. As the magnitude of (gGM)r continues 
to fall below that of {gGE)r, it reaches a minimum at t ~ 0.3 to in the galaxy evolution phase, and then once again climbs 
above that of ((?GE)r at t ~ 0.9 to as galaxies continue to evolve, reaching a value of ~ 1.4 x 10~^ cms"'^, greater than that of 
(pge))- (~ 1.7 X 10~* cms^^), at t = to, indicating the possibility for the Universe to be in an accelerating phase for ^^0.5 
[compare Figs. 2(a) through 2(e)]. 



5.2 The GM and the GE Accelerations at the Hubble Radius: From Time of Planck Scale to Inflation to 
the Present 

Upon substitution of equation (|54p into equation (I48p . setting r ^ 0, 6 = ^ , and using model parameters defined in Section [331 
(eauation l43l a = exp[— 115(t/to)], k = 71a, Ht ~ 1/t), we can derive an approximate analytical expression for the GM cosmic 
acceleration that appears to be valid at early times in the Universe as well as later times when the scale factor of equation (|43[) 
is normalized at the present epoch [i.e., when fi of equations 1)50^ to 1)52^ equals 1]. Thus, we find that 

/ N 1 o in7 *o / 57.5t \ cm 

(<7GM)r ~ 1.8 X 10 ^ exp I I — , (55) 

where n is defined in equation (|43p . For t — to (i.e., H — Ho) equation (|55|) gives the same value as that calculated using the 
exact analytical expression for the GM acceleration (eg nation 148 II . measured at an arbitrary comoving present epoch observer: 
(<7GM)r- ~ 4 X 10^* cms~^, for Ho = 71 kms~^ Mpc~^, by letting z — > [or r — > 0; compare equation 1521 and Fig. 2(f)]. 
This appears to validate our use of equation (|55|l at earlier times for which r ~ and t 7^ to. 

Displayed in Fig. 4 are the evolutions of the magnitudes of the GM and the GE accelerations ((?GM)r and ((?GE)r, 
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respectively. Plotted in Fig. 4(a) are the magnitudes of {gGM)r (equation I55[) and for!F,)r (equation I49|) . from the Planck time 
{tp ~ 5.3 X 10""" s, with length r = = ~ 1.6 x 10"^^ cm) up to the Hubble time of t ~ 4 x 10^ s (~ 138 yr = 10"** to)- 
Note, for any given epoch, the magnitudes of {gayijr and ((?GE)r are evaluated at the Hubble radius, i.e., the limit of the 
causally-connected observable Universe for any arbitrary observer. If we assume that n = | , being consistent with a Lemaitre 
early cosmology, then from the Planck time up to the time of inflation (~ 10~^^ s), {gGM)r goes from ~ 10^* down to ~ 10^^, 
respectively, factors larger than |(<?GE)r|. The net acceleration of the expansion might be suflicient to inflate the Universe to 
outside of a causally connected region: given by the Hubble radius (i.e., vh = ct = cH~^), consistence with what we assume to 
have occurred in our standard FLRW cosmological with inflation. The basic requirement is that the scale factor must increase 
by a factor of at least 10'^^ in a very short time. This could imply that perhaps what we call inflation may have been, or at 
least in part, caused by {gaM)r, and thus is related to frame dragging in some degree. The above statement, however, would 
have to be investigated. Notice, even up to age f ~ 1 s, the proposed limit for inflation to have occurred to be consistent 
with nucleosynthesis (Liddle 2003), ((?GM)r is still ~ 10^^ orders of magnitude larger than |((;GE)r|. Since, however, we want 
this present model to be consistent with the standard cosmological model, and with the assumption that {gGM)r is related to 
spacetime torsion, we must keep in mind that the negative terms in the last acceleration component on the right-hand side 
of equation ((5| could very well come into play to keep the expansion of the Universe at a rate consistent the standard model. 
Importantly, ai t — 138 yr {t ~ 4.4 x 10® s), ((?GM)r has fallen signiflcantly: (gGM)r ~ 41.5 cms~^ and ]((7GE)r| ~ 1.7 cms~^, 
for n = |, implying radiation dominance, n' — 1 (compare equations 1431 and 144 |l . 

In Fig. 4(b), we continue the calculations of cosmic gravitational accelerations of the GE (equation I48|) and GM (equa- 
tion l49p radial components over time from t ^ 138 yr up to the present, still evaluated at the Hubble radius; rnit) = cH^^ , 
the causally connected region surrounding an arbitrary comoving observer who is at the centre (r = 0) in the metric of equa- 
tion ([7|, with this centre being any point in spacetime. The proper physical distance r{t) (equation I37p . measured according 
to a standard clock at rest with this arbitrary comoving observer, for r{t — to)H = cHq^ , is equivalent to redshift « ^ 30 
according to the relativistic Hubble law: 



V 



(1+^) 



-no 



for Ho = 71 kms ^Mpc ^, where z can be solved for: 

/ l + rHo/c ] 
[l-rHo/cJ 

Note, r = rn for z ^ 30, according to equation (|56p . appears to simply suggests that r 5C rn is always the observable causally 
connected Universe for any z ^ 30, which is relativistically true. Equation (|56|) gives the spatial separation at a common time 
from an arbitrary comoving observer when the light was emitted from say a distant galaxy at r ^ r^f. 

In other words - making a brief deviation to clarify this cosmological distance given by equation ((7|) - the relativistic 
Hubble law, gives the distance measured by causally connected observers at a common time t. Simultaneity, means setting 
dt = in the metric of equation ((7| (or in the standard Friedmann- Robertson- Walker metric), which implies choosing the 
local frame of this freely falling comoving observer. The proper distance between spacetime events is then given, along the 
line-of-sight (i.e., in the radial direction, with dcp — d9 — 0), by 



\/\9rr\dr, (57) 







where x is the comoving coordinate distance (associated with the proper physical distance), discussed in Section [3.41 and it 
can be shown from the metric component, grr, of equation ([8]), upon integration, that 

^proper = R{t)X = 



which is just equation (|37p . The proper distance measured by a comoving observer can be understood as a spacelike separation 
using a hypothetical ruler to measure the separation at the time of emission, from say a distant galaxy, as opposed to a lightlike 
(null) separation using light-travel time in an expanding universe to measure the separation. Both r{t) and x &re spacelike 
separations between events: this means that they are imaginary (oc i = 1) and cannot lie on the world line of any body or 
particle. From the proper distance r{t) of equation H37[). its derivative with respect to time, and equation ([6]), the Hubble law 
can be derived exactly. The proper distance might be called the dynamical distance. The recession velocity for the proper (or 
so-called dynamical) cosmological distance is always ^ c, at r ^ vh, i.e., causally connected regions. Now, since the Hubble 
law is derived exactly from the metric proper radial distance, and since the dynamics are what we are concerned with in this 
present manuscript, we appropriately use the proper distance in these calculations. 

Finally, before proceeding, for the proper distance r, if we consider the ratio of the approximated distance (equation [52} 
to the more accurate relativistically corrected distance (equation [56}, we obtain for z — 0.03, 0.5, and 1 the ratios ~ 
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1, 1.3, and 1.7, respectively. In these present model calculations, to give an explanation for the observed so-called dark 
energy, we need only consider the evolution of the GM (equation [48]) and GE (equation [49]) radial accelerations over spacetime 
points for z < 1, since equation (|48|) goes to zero, as measured by a present epoch comoving observer, for z « 0.7 using 
equation (|52[1 . the nonrelativistic Hubble law; and, it goes to zero for z ~ 1.5 using equation (|56p . the relativistic Hubble 
law [compare Figs. 2 and 4(b)]. Both these z values where (pGM)r — > are qualitatively consistent with observations. Thus, 
the results of this present model, giving an explanation for so-called dark energy, do not change qualitatively, and these 
results change only slightly quantitatively, when the approximate distance f equation I52|) . as apposed to the relativistically 
corrected distance (equation [56|, is used for z ^ 0.68. Importantly, the accelerated expansion appears to set in, theoretically, 
according to this model, at 2: ~ 0.5, consistent with observations, when using either the nonrelativistic or relativistic Hubble 
law to determine r, the physical distance. Therefore, the above, and the fact that most large-scale galaxy surveys use the 
nonrelativistic Hubble law, appear to justify our use of equation H52|) in this present manuscript. Nevertheless, some relevant 
results using the more accurate relativistic Hubble law will be discussed below. 

Proceeding, the early evolution {t ^ 10~*to) of the GM acceleration as measured by an arbitrary comoving observer at 
a specific time is expected not be qualitatively, and perhaps not much quantitatively, different from that seen in Fig. 4(a), 
where {gGM)r in these times is given by equation (|55p for sufficiently small r. Thus, in these early times {gGM)r decreases 
with increasing time or decreasing angular velocity, as would be expected (see also Fig. 3), independently of r, recalling that 
r = r{t). But the evolution of {gGM)r (equation I48|) from t — 10~* to {~ 138 yr) up to the present epoch will be quite different 
from that of Fig. 4(b) {z ~ 30), for a sufficiently smaller value of z (associated with an r), as can be seen in Fig. 2, wherein 
cosmic acceleration of the expansion appears to begin at z ~ 0.5, as measured by a present epoch arbitrary comoving observer; 
this can be seen in Fig. 2(c). We discuss this further in the next paragraph. 

Displayed in Fig. 4(b) are the magnitudes of the accelerations, ((?GM)r and ((7GE)r, between the interval 138 yr ^ t ^ to, 
measured by an arbitrary comoving observer, at the Hubble distance r(t) = ru = ct, i.e., at 2: ~ 30 according to equation (|56p , 
as this spacetime coordinate distance (or point) evolves over time, with to indicating the present epoch Hubble radius (causally- 
connected observable Universe). In this case, {gGM)r falls below the magnitude |((;GE)r| twice, at t ~ 2 x 10~^ to and t ~ 
5 X 10~^to, because of the so-called artifact in changing from radiation dominance to matter dominance occurring at the 
equilibrium time t = teq ~ 4 x 10~^ to. This is where n — > | at the radiation-matter equilibrium time t = tcq — 1.7 x 10^^ s 
(= 4 X 10"** to), for fimat =i 0.3 (Liddle 2003); compare Fig 4(b). Then, although not shown on Fig. 4(b), at t ~ 0.6 to 
(~ 8 X 10^ yr), ((?GM)r falls to zero, while |((/GE)r| ~ 5.8 x 10~*cms~^. After this, {gGM)r remains zero up to the present 
epoch (to) at which |((7GE)r! ~ 3.5 x 10^*cms^^, consistent with a decelerating universe, as measured by a present epoch 
arbitrary observer. But, recall, as z becomes less than one, at t = to, as measured by a present epoch observer, {gGM)r becomes 
larger than ((/GE)r for z ^ 0.5, indicating perhaps the Universe enters into an accelerating phase [compare Figs. 2(a) through 
2(e)], consistent with recent observations. Note, as mentioned earlier, at 2 = r = and t = to, indicating at the present epoch 
arbitrary comoving observer, ((/GM)r ~ 4.9 x 10^^ cms~^ [compare Fig. 2(f)]. 

After a quantitative comparison of the calculated results, it appears that [gGM)r is independent of r for small r, i.e., 
(<7GM)r has the same value for any value of r(t) as measured by an arbitrary comoving observer at a specific epoch up to 
some time t = tcrit. For example, notice in Figs. 2 and 4(b), over the range of values of z measured by a present epoch 
arbitrary observer, the value of {gGyijr — 41.5 cms^^ at t = 10~* to is the same, although the values of r(t) are different. 
Analytically, this is because the exponent of the exponential term in equation (|48|) goes to zero, meaning e'' — 1, for small r; 
and because (m/c)rsin6 is <^ 1, where m is given by equation Q; thence ((Jgm))- appears independent of r (or 2); compare 
equations (|48p and 1)55^ . Now, comparing different values of 2 in the range 4 x 10~® 5; 2 ^ 30, with {gGM)r evolving over 
time, as in Figs. 2 and 4, with the scale factor being normalized at the present epoch, it appears that {gGM)r begins showing 
dependence on r(t) at the critical time tcrit — 1-2 x 10* yr (= 9 x 10~^to) after the Big Bang, and then showing significant 
dependence at times t ^ 2.8 x 10* yr (= 0.02 to). This means that {gGM)r begins to change significantly over time, wherein 
it starts to increase. Before discussing below the importance of this observation, we readily see that this validates our use 
of equation (|55p to evaluate {gGM)r in the early Universe: it is independent of r for small r and it matches or converges to 
the value of {gGM)r — 41.5 cms~'^, evaluated using the exact analytical expression (equation [48]), at the so-called interface: 
t = 10"** to = 138 yr ~ 4.4 x 10^ s (compare Figs. 2 and 4). 

Continuing, this independence of ((?GM)r for small r, at early times, as the Universe evolves over time, and then becoming 
significantly dependent on r (or 2), as measured by an arbitrary comoving observer, for t ^ 0.02 to, can be seen in Figs. 2 
and 4. This behavior can be further understood by looking at the resulting equation (|55[) for very small r at early times 
and for very small r at later times. Equation 1)55^ shows what happens for small spacetimes, when the exponential goes to 
one in equation (|48[) . as mentioned above, and ((Jgm))- decreases somewhat linearly with increasing time [compare Fig. 4(a)]. 
But as time continues to increase, the exponential in equation (|55p begins to increase, and {gGM)r reaches a finite value at 
t = to, 2 = 0, like that seen in Fig. 2(f), which was evaluated, along with the other figures of Fig. 2, using the exact analytical 
expression for ((/GM)r (equation |48|. So, a present-day arbitrary comoving observer will measure the value of ((?GM)r to be 
different at different 2 (or r) values; and, importantly, its values at 2 ^ 0.5 dominate over |(3ge)i1, contributing possibly to 
the recently observed acceleration of the cosmic expansion. 
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Thus, looking at Fig. 2 and what happens sA, 0.02 to enable us to see how the GM acceleration {gGM)r might once 
again dominate over the GE deceleration {gG¥^)r, if it were an intrinsic part of the equation of motion of the expansion 
(or cosmic scale factor) in the expanding and rotating Universe of equation ([5]) to produce the observed present-day cosmic 
acceleration. This we suppose in Sections 15.51 and 15.61 where we will consider, as mentioned earlier, how the GM acceleration, 
and thus inertial spacetime frame dragging, might be related to the torsion term in the Godel-Obukhov equation of motion 
(equation [5j . 

Further, and for completion, from using the relativistic Hubble law (equation I56|l to determine r, the physical proper 
distance, in equation (148 p . the following appears to occur. For any 10, (gGM)r — > (or negligibly small) at around 
t = 0.6 io — 8.3 X 10^ yr after the Big Bang, i.e., about 5.5 x 10^ years ago, and remains zero up to the present epoch 
to = 13.8 X 10^ yr (compare Fig. 4). There is no turning up of the curve, like that of the curves found in Fig. 2, of (pGM)r versus 
cosmic time, for these large values of z. The case of Fig. 4(b), which is evaluated at the Hubble radius rn, equivalent to 2: ~ 30 
according to the relativistic Hubble law, is an example of this. Then from around z = ((?GM)r gradually changes from zero to 
nonzero values over time and decreasing redshifts, as measured by an arbitrary comoving observer, consistent with observations 
that indicate the coming presence of so-called dark energy. At 2: ~ 1 is where ((7GM)r ~ 2.4 x 10"^^ cms^^ first becomes 
nonzero at t = to, as measured by a present epoch observer, but {goyCir is still much less than |(gGE)r| ~ 2.1 x 10~* cms~^ 
at this redshift. Then at around z — 0.5, as mentioned above, ((?GM)r becomes greater than |((;GE)r|, indicating the Universe 
goes into an accelerating expansion phase [compare Fig. 2(c)]. Further details of these results using the relativistic Hubble 
law, at these large redshifts, are presented elsewhere (Williams 2011, in preparation). 

So, in summary, the above appears to suggest that the GM field strength was once very large in the past as shown in 
Fig. 4(a); became negligibly small (~ 0) around t — 0.6 to, about 5.5 x 10^ years ago, for z^ 10, as shown somewhat in 
Fig. 4(b); and, then, from z ~ 6 the GM field gradually increased over time, and decreasing z, from zero to a nonzero value 
at 2 ~ 1, as measured by a present epoch observer at t = to. This is the first indication of a turning up of the curve (compare 
Fig. 2 and Fig. 4). Finally, at around 2 = 0.5, the magnitude of (pGM)r became greater than that of ((7GE)r and the Universe 
entered into an accelerating phase, thus, consistent with recent observations. 



5.3 The GM Acceleration, Cosmic Rotation, and the Present epoch Observer 

Figure 5 displays what happens to equation (|48p in a general sense when 2 is held constant ai z — 0.5 [Fig. 5(a)] or at 2 = 
[i.e., r = 0; Fig. 5(b)] and 5 is allowed to vary, where 

(a;,.ot)o = 5 Ho (58) 

is the present-day angular velocity (see eauations l40ll41l and 142 p . with Ho — 71 kms~^ Mpc~^ = 2.3 x 10~^* s~^. In Figure 5(a) 
6 was allowed to vary between 74 ^ 5 ^ 1"^**. We see that ((?GM)r — > exponentially (i.e., quickly), left of the maximum, at 
any value (ajrot)o ^ IOHq: falling from a maximum of ~ 5 x 10^^ cms~^ sd, 5 — 0.1; and ((;gm)j- — > linearly (i.e., slowly), 
right of the maximum, as (ajrot)o — > at S = 10~'^*. Note, (a;rot)o = 2tvHo — 1.5 x lO^^'^ s^^ is the value used in the 
calculations of this present manuscript (see equation I54|) . The reason for this behavior is that, for a present epoch observer, 
with R{t = to) — 1, at large S, the exponential expression in the denominator of equation (|48p dominates, causing ((?GM)r 
to go exponentially to zero; and at small 5, the exponential expression goes to one and (aJrot)o in the numerator dominates 
(compare equation |4]) , thus, causing {gaM)r to go to zero in a linear-like fashion. This behavior will be important when we 
next analyze the equation of motion describing the expansion and deceleration of the Godel-Obukhov spacetime cosmology. 
We shall see in the following section that this behavior is similar to what one would expected of the third term in equation ([5]), 
again with R{to) = 1. Figure 5(b), with 5 allowed to vary between 74 ^ 5 ^ 1~^ at r = 2 = 0, shows, again in a general 
sense, what the strength of {gaM)r would be as measured by a present epoch observer (t = to). The reason for the linear-like 
decline in {gGM)r is because with r = the exponential expression in the denominator of equation (|48p equals one, resulting 
in {gGM)r decreasing as (a;rot)o in the numerator decreases. 

Although the above analysis presented in Fig. 5 may not have any real physical significance, save limiting the proportion- 
ality constant relating Ho and (ajrot)o (equation [58J , it shows the general behavior of the dependence of {gGM)r on (a;rot)o at 
2 = 0.5 and 2 = 0, as measured by a present epoch observer. 



5.4 Comparing the Standard FLRW Spacetime with the Godel-Obukhov Spacetime 

In this section we will analyze the equations of motion of the scale factor that contains terms that control spacetime acceleration 
and deceleration of the Universe over time: equations (|32p and ([5|, for the standard FLRW and the Godel-Obukhov spacetimes, 
respectively. Substitution of equation ((33]) into the second term on the right-hand side of equation (|32p . and comparing 
equation (|36p . allows us to identify this second term as the GE acceleration that decelerates the Universe, particularly when 



Could Dark Energy be a Manifestation of Gravity? 21 



E 




1e-20 1e-25 1e-30 1e-35 1e-40 
omega (1/s) 




1e-16 1e-17 1e-18 1e-19 1e-20 
omega (l/s) 



Figure 5. The GM cosmic acceleration and rotation: (gGM)r vs. uiiot- (a) At redstiift z = 0.5, scale factor R{t) - 
5 in equation II58I I decreases along the curve, as LOrot decreases (see text), (b) At redshift z = 0, scale factor R(t) 
for comparison, using 5 = 2n (see equation I54II . as it is done throughout these calculations, at t = to, with Hq 
(aJrot)o ~ 1.5 X 10-" s-i. 



1, at t = ig. Note, 
1, at t = to. Note, 
: 71 kms-^Mpc-l, 



I = -qH', (59) 
or 

1 = gpLRW = ^^JJ2- (60) 

Equation (|60l) is just equation U]), where equation ((6]) has been used. Notice, we have defined q = gpLpivv distinguish between 
deceleration parameters in the two spacetimes we are considering: the FLRW and the Godel-Obukhov, where q = q^^ in the 
Godel-Obukhov spacetime (indicated by the subscript 'GO'). We will see below how the two may correlate. So, we find in 
general from equations (|59p and (|60|) that 

9flrw = =^ coasting, 
9flrw > deceleration, 

'J'flrw < acceleration, (61) 
again, for A = p = 0. 

On the other hand, in the Godel-Obukhov spacetime we can identify the GE deceleration as the first term on the right- 
hand side of equation ((5]), and exactly that of the standard FLRW model if we assume q = 1 (see equation [59]), as mentioned in 
Section [S] Moreover, it seems reasonable to identify the fourth term on the right-hand side of equation ([5]), involving the mass 
density pt, as that describing the initial inertial expansion over time. We define this acceleration term, which we associated 
with the initial cosmic expansion by 




k 



k 



qH\ (62) 



where equation (|33p has been used (see Section [3.41 for validation of its use). This term being proportional to and, thus, 
the density, will be very large in the very early Universe. Now using equation (|37[l and its derivatives, we can express this 
inertial spacetime cosmic expansion in terms of the physical separation distance r: 



k + o 



r = {ai)r = 2 qH^r, (63) 

where r is given by equation (|37p . So, we have identified equation (|63[) as the acceleration due to expansion. If we use the 
expressions for o and k given in equations ((38} and (|39p . respectively, with values given in Section [3.51 for constants ci and 
C2, equation (|63p reduces to 



{ai)r~2qH'r. (64) 

This is consistent with the of order expression given near the end of Section [2] which was derived from dimensional analysis. 

At this point we will used the analogy of the standard FLRW cosmology, where we set K = p = K (curvature) = 
in equation p2p to define the deceleration parameter (Qflrw) ''■^'d its relation to the mass density pt; here we will set 
i^rot = B = p = in equation ([5]) to define the deceleration parameter (ggo ) and its relation to the mass density. Equation ([5]) 
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reduces to 

R „2 , 8nG f k + a , 



(compare with equation 1321 using equation [33]). Using the model parameters of Section [3.51 as was done in equation (f64l 
substituting in the mass density now given by 



and dividing through by —H^, equation (|65|) yields 



(66) 



.oo4(^^ + lj; (67) 
or 

(compare equation I60p. So by equation (|68p we now have a relationship between the deceleration parameters of the standard 
FLRW and Godel-Obukhov cosmologies, which appears to be consistent with observations, as can be seen in the following: 
For 

g^o ~ i, gpLRw = 0^ coasting; 

< i, gpLnw > ^ deceleration; 
(Jqo > i, gpLRw < ^ acceleration; (69) 
compare equation (|6ip . 

Next, looking at the second and third terms on the right-hand side of equation ((5]), it appears that these are related to 
the fictitious-like forces associated with cosmic rotation, and at least one to general relativistic inertial frame dragging. These 
fictitious forces, analogous to Newtonian centrifugal and Coriolis forces, appear in the equation of motion of an object in a 
rotating frame. It is called a fictitious force because it does not appear when the motion is expressed in a Newtonian inertial 
frame of reference (i.e., a frame that is not rotating nor dragged into rotation). The second term can be easily identified as a 
centrifugal-like acceleration, associated with the vorticity (i.e., the rotation), that decreases over time, being proportional to 
LJrot, and largest in the early Universe. This second term appears to be associated with the initial cosmic rotational energy, 
similar to the fourth term on the right-hand side of equation ([5]) in which we identified above (equation I62|) as that being 
associated with the initial inertial cosmic expansion energy. The third term, as mentioned earlier in Section 13.11 is related 
to the torsion, coupled with spin (or rotation) and curvature, of spacetime, and appears to be directly related to general 
relativistic inertial spacetime frame dragging, where we would expect it to be some sort of Coriolis-like force. This third term 
behaves similar to that seen in Fig. 5(a), i.e., it goes to zero for large Wrot, which would be at early times; increases as ojrot 
decreases over time; and, then, at later times, decreases as, perhaps, the magnitude of the cosmic magnetic field B (Obukhov 
2000) decreases over time. Note, we shall consider this B and its relation to Wrot further in the following section. Now, since 
in the case of a moving test particle in the gravitational potential well of a compact rotating object, the GM force exerted 
on the particle has been related to a Coriolis-like force, one would expect the GM general relativistic so-called fictitious force 
(exerted on a moving 'test galaxy') here, and thus the third term on the right-hand side of equation ([5]), to be related to a 
Coriolis-like force. So, supposing this to be the case, based on 

(i) the evolution of {gGM)r presented in Fig. 2 and its consistency with recent observations of cosmic accelerated expansion 
as discussed in Sections 15.11 and 15.21 

(ii) that (gGM)r is a Coriolis-like force derived using the Godel-Obukhov metric (equation[7|, and 

(iii) support by the behavior of the third term: being somewhat similar to that of Fig. 5(a), 

it seems reasonable that the acceleration produced by this third term should be at least approximately equal to that of the GM 
force given by equation (|48|) . But, in order to set these accelerations equal, and solve for the cosmological coupling constants, 
Ai and A3, whose derivations are otherwise very complex (see Obukhov 2000), we must first express equation H48p in the form 
of a component of the cosmic scale factor equation of motion (equation [SJ. This is done in the following section. 

5.5 The GM Field, Spin, and the Electromagnetic Field 

In this section we will analyze the terms in equation ((5| involving the magnetic field, B, and the meaning of this electro- 
dynamic field. As a vector representation this field is the F12 — B^ component of the electromagnetic field tensor due to 
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the electrodynamic characteristics of the spacetime matter. Here we consider the Universe to be a spinning Einstein-Cartan 
(source of spacetime curvature) fluid of charge density rotating about the axis, i.e., ojrot = tt'rotez, like that stated in 
Section [3] (Obukhov 2000; Obukhov & Korotky 1987). We assume that torsion is generated by the spin tensor of such a fluid. 
Importantly, and somewhat surprisingly, it appears that B is related to the GM field, as we shall see below. 

Using equation (|37|) . its derivatives, and dividing equation (|48l) through by r, we have the GM acceleration expressed in 
referring to the above section, in units of equation ((5]), i.e., s~^, expressed as an acceleration of the cosmic 



desired form 
scale factor: 



GM 



ka{t) 



[k + a{t)] 



1/2 



{m/c)r sin 9 + 1 
ii:(t)sin6le(™/'=)''='" 



a{t) + 2k \ LJrotC 



2[k + a{t)] 



(70) 



where, m, being proportional to ojrot, is again given by equation Q. Notice that the terms in the above equation are consistent 
with the GM force producing a Coriolis-like acceleration. 

Note, in these calculations (as mentioned in Section |4] in reference to equation [48]) , we set 6 — ^, measured locally by 
an arbitrary comoving observer. This appears to be consistent with the quasi-Cosmological Principle in a rotating universe, 
where, for this local observer at the centre of the metric (equation[7|, we propose four-dimension spacetime to be described (or 
at least understood) by two-dimensional spacetime spherical axisymmetrical concentric hypersurfaces of spacelike separations. 
We assume each surface, surrounding the local observer, to be embedded with the constituents of the Universe at any given 
time, and rotating at the corresponding angular velocity, cjrot = iJTot{t) = (cJrot)t (compare equations 1541 and I53 [) . Each 
surface of infinitesimal thickness represents a specific epoch of time, with time increasing (locally) toward the centre of these 
concentric surfaces of constant time t as well as constant tJrot, for observers whose world lines coincide at any point on a 
particular surface. An arbitrary comoving observer of the metric of equation ((7]), with ^ = f , would be at the centre of a 
cross section of such axisymmetrical and homogeneous surfaces of constant time. It might be easier to visualize this local 
arbitrary observer to be at the global 'centre' of the Universe (although not physically nor a necessity) such that the local 
coordinates coincide with the global system. So this arbitrary comoving observer's line-of-sight at any given time and location 
will be from the centre of a disc plane of infinitesimal thickness, intersecting surrounding concentric hypersurfaces, whose 
radial distances from the observer are measured by r{t), the proper or physical distance (see equations [37] and [57]). Such 
a geometrical configuration for a local arbitrary comoving observer with 9 = ^ appears to be unique to local observers in 
a rotating universe, and appears to be consistent with the quasi-Cosmological Principle of homogeneity and isotropy in the 
CMB radiation, but now including spatial isotropy in the direction of r (for axial symmetry), i.e., pseudo-spatial isotropy in r 
resembling that for the nonrotating FLRW spherical symmetric case. This proposed geometry is for a present-day comoving 
arbitrary observer, where the scale factor R{t) (equation [43]) is normalized at the present time t — to, at r = 0, corresponding 
to redshift z = 0, with {gGM)r ~ 4 x 10^* cms~^ at to, as can be seen in Fig. 2(f). As r, the proper or physical distance, 
increases, this observer is looking back in spacetime to an earlier time, indicated by concentric spacetime hypersurfaces. This 
is the arbitrary comoving observer that measures the parameters displayed in Figs. 2 through 6, described in Sections [4] and[5] 

Now with the note above in mind, we set the third term on the right-hand side of equation ([5]) that is associated with 



torsion, approximately equal to that of equation (|70[1 above, associated with frame dragging, and solve for 4A3 — A^, the 



coupling constants of torsion, curvature, and spin, 
suggested to be the case in recent observations: 



with 4A3 2> Ai to prevent cosmological collapse (Obukhov 2000), as 



4A3 - A? 
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k + a J \ [fc + cr(t)]3 
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(m/c)r sin 6^ + 1 
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1 - 



a{t) + 2k \ R'uj^otC 



2[k + a{t)] 



B^r 



(71) 



where in cgs units, Ai and A3 have units of cmg~^, and B has units of gauss {— g^'^^cm~^''^s~^). The torsion, as related to 
torque and rotation, appears to be caused by the inertial spacetime expansion as it attempts to drag (or torque) the cosmic 
spacetime matter out of rotation (as described in Section [2]). I refer to the cosmic expansion frame as an inertial (or 'flat') 
spacetime frame in a general relativistic dynamical sense because it has inertial force properties as well as inertial motion 
properties, appearing to be, in straight lines away from an arbitrary comoving observer (i.e., along the observer's line-of-sight) 
as expected of an inertially flat expanding spacetime frame. In other words, in a general relativistic sense, the inertial frame 
of linear expansion is torquing or dragging the non-inertial frame of rotation, of cosmic spacetime mass density pt, out of 
rotation, in a sense, producing the GM force (per unit of moving mass) that accelerates the cosmic expansion. So, the torsion 
it appears is caused by the expanding frame's inertia to rotation, producing the negative frame dragging angular velocity (see 
equation II 1|) . This behavior is analogous to how a sufficiently large mass density can warp (or curve) spacetime causing the 
inertial cosmic expansion to decelerate, clearly seen in the standard FLRW cosmological model (equation [32]), thus producing 
the so-called GE force (equation [31] compare also equation [5]). 

Moreover, in general relativity, the effect that frame dragging has on moving objects in spacetime is described by the GM 
force field. Similar to that experienced by moving objects (i.e., test particles) in the gravitational potential well of a rotating 
black hole [see fig. 2 in Williams (2002)], where the frame dragging angular velocity is in the positive azimuthal direction and 
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produces a positive radial force, this cosmological spacetime frame dragging is, however, in the negative azimuthal direction 
but produces a positive radial force as well. The difference in the sign of the frame dragging velocity appears to be due to the 
involvement of both the translational inertial and rotational non-inertial frames in an expanding and rotating universe, acting 
on freely falling or comoving observers, such as galaxies; whereas, with the black hole only the rotational frame produced 
by the black hole is involved, acting on freely falling local inertial frame observers, such as particles of plasma, dragging 
these observers into rotation in the direction that the black hole is rotating. Now, in general, both the translational cosmic 
spacetime inertial and cosmic rotational non-inertial frames are undergoing a form of acceleration: a change in the velocity for 
the translational or expansion frame and a change in both direction and velocity for the rotational frame. Then according to 
Einstein's Equivalence Principle of Gravitation and Inertia both frames give rise to gravitational accelerations (attractive or 
repulsive) such as those found in equations (O and (|32p . Finally, we point out that angular momentum of the cosmic rotation 
appears to be transferred to linear momentum of the inertial expansion, and conserved, as the expanding spacetime frame 
appears to account for the angular momentum lost to torque, thus, maintaining a close (or harmonic) relationship between 
Ht, the expansion rate, and {LOrot)t, the rotation rate (equation I54p. as suggested by observation (equation I40p. 

Next, the derivation of equation (|7ip suggests that there might be a relationship between the GM and the electromagnetic 
fields; we will look further at this somewhat surprising discovery below. The magnetic field given by Obukhov (2000) is from 
application of an ideal fiuid plasma of spin and energy-momentum, to a cosmological model with rotation (shear- free) and 
expansion, in the framework of the Einstein-Cartan theory of gravity. Note, the Einstein-Cartan theory of gravity is just 
Einstein's theory including rotation (and its effect on spacetime), a natural extension to describe a universe with expansion 
and rotation. This cosmic magnetic field, seen in equations ^ and (|71|l , and mentioned above, is the F12 = Bz component of 
the electromagnetic field tensor describing the electrodynamic characteristics of matter in the Universe. It appears that this 
field can be expressed over time by (Obukhov 2000) 

Bt = [-2R{t)uj,^trt]^, (72) 

where rt [= r(t) < 0] is the spin density, in units of gcm~^ s~^\ again, R(t) and ojrot = (tt^rot)t are given by equations (|43|l 
and (|54p . respectively. The spin density by definition is the angular momentum per unit volume. It then seems reasonable to 
express the spin density as 

Tt ~ —ptUJrotr'^, (73) 

where, also, r = r{t). This means that B could be a primordial cosmic global magnetic field intrinsic to the rotating and 
expanding spacetime matter in the Universe. Since B of equation (|72p depends on the spin density, rt, and this spin density 
can be expressed in terms of the average mass density pt by equation (|73|) . with great surprise and excitement, equations (|72|) 
and (|73p appear to show or at least lead us to the long sought relationship, by Einstein (as well as by the author after 
being inspired by the work of Einstein), between the electromagnetic field (with the electric field canceling assuming equal 
number of positive and negative charges) and the gravitational field, namely, between the magnetic field and the GM field. 
It appears that they can be related by the spin density. Further details of this relationship and of its importance in other 
astrophysical phenomena are beyond the scope of this paper and will be discussed elsewhere in a forthcoming paper (Williams, 
in preparation). But important here is that these findings, in the context of the model considered in this paper, reveal how 
the GM field might be associated with the magnetic field through the spin density rt, which can be expressed in terms of the 
average mass density pt- 

In the following we shall compare values given by equations (|72|) and (|73|) in the early and late Universe with other 
theoretical estimated values, as well as with observations. First, for comparison, we used the author's derived spin density 
Tt of equation (|73|) and the cosmological parameters used [qo — 0.01 and (airot)o = 0.1-ffo] by Obukhov (2000), with Ho ~ 
71 kms~^Mpc~^, to see if the values of the author and Obukhov (2000) agree for spin density, ro, at the present epoch. The 
density from equation (I34p is calculated to give pc — 1.89 x 10""^^ gcm~^; then, evaluation of equation (|73|l at the Hubble 
radius (th — 1.3 x 10^* cm) yields ro ~ —7.4 x 10^ gcm~^s~^. The magnitude of this value is smaller than the magnitude 
of Obukhov's (2000) estimated value (ro ~ —5 x 10* gcm'^s"'^) when using the same go and (ajrot)o as used by Obukhov 
(2000). However, using the same qo — 0.01, but letting (a;rot)o = 2ttHo (equation [54]), as used by the author, equation (|73[) 
then yields ro ~ —4.7 x 10* gcm~^s~^. As we can see, this value is approximately equal to the average value estimated by 
Obukhov (2000), given above. Since Obukhov (2000) does not give the derived analytical expression used to calculate ro, it 
cannot be commented on here as to why our values agree for the same qo, but with a different relationship between Ho and 
(i^rot)o, except that perhaps this is the reason Obukhov (2000) calculates the modern-day magnetic field strength to be ~ 10^ 
times larger than the established upper limit from astrophysical observations. Nevertheless, from the above, it seems safe to 
conclude that the deduced analytical expression for rt (equation I73[l is at least of order true. 

Next, the cosmic magnetic field of equation (I72p can be expressed as 

Bt '-^ [2R{t)u;l,ptr^]i , (74) 
where we have used equation (|73p . It appears that Bt above is a frozen-in cosmic primordial magnetic field modified only 
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through the expansion process, such that the magnetic flux is conserved (Asseo & Sol 1987). Then Bt would be that measured 
by a comoving arbitrary observer at the centre of the metric of equation ([7]). As the physical distance r from this observer is 
increased, the observer is looking back in time, as usual, because of the finite speed of light. So with this in mind one would 
expect Bt to be larger at large r, as measured by a present-day observer (i.e., where the scale factor of equation|331is normalized 
to one), and smaller as r gets smaller, until it reaches the value measured locally, between galaxies. This Bt = Bq would be 
the present-day value of the cosmic magnetic field that has been somewhat dissipated, due to the spacetime expansion of the 
Universe. Below we will calculate values measured by this arbitrary comoving observer. 

At the present epoch, as measured by this comoving observer, with q = i, n = | using equation (|43|) . normalized at 
t = to, po = Pc{to), for Ho = 71 kms"^Mpc"^ = 2.3 x 10"^* s"\ luq ^ 1.5 x 10"^'^ s"^ (equation [54)) . to = H^^, and 
r — rn = cH^^ , we find that _Bo ~ 8 x 10~* gauss at the Hubble radius (~ 1.3 x 10^* cm). This can be seen in Fig. 6 where 
the magnetic field of equation (|74[) is evolved over time from when the Universe was 138 years-old to the present. Figures 6(a) 
and 6(b) display the evolution of the field strength Bt measured at the Hubble radius {rn, at z ~ 30 according to eauation l56|l 
and at z = 0.5, respectively, by an arbitrary comoving observer. These figures are to be compared with those of Figs. 4(b) 
and 2(c), respectively, which plot the GM acceleration over time. 

For the very early Universe, this comoving observer measures the following as spacetime evolves. At the Planck scale 
(indicated by subscript P), t = tp = ~ 5.4 x lO"*"* s, with pp ~ 1.8 x 10^^ gcm"^ (equation [M]), Bp ~ 

1.8 X 10'^^ gauss, according to equation (|72[) . at the Hubble radius (r = rp = ctp ~ 1.6 x 10"'^'' cm or the so-called Planck 
length), where h is the Planck constant. The best way to express Bt, it appears, at least during inflation, and to see clearly 
how the magnetic fleld falls off over time, is to use 
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then equation (|74|) can be expressed as 

(76) 



Bt = 



2ajrotPor^ 



Now, we assume that inflation occurs at the characteristic times between 10~^® s ^ finfl ^ 10"'^* s, where during inflation 
Ht ~ constant, since At ~ 9.9 x 10~^^ s <C 1 (see Section [3.5p . The Hubble radii (r — rinfj ~ ctinfl) corresponding to the 
time interval above are 3 X lO-^*^ cm < Hnfl ^ 3 X 10"^* cm. From equation (|74p or equation (|76|) we calculate the magnetic 
fleld at the beginning of inflation to be Bt ~ 4.7 x 10^^ gauss. Since we assumed that Ht is approximately constant during 
inflation, it seems reasonable to assume that Urot (~ 2ttH) is also approximately constant. We found in Section [3.51 that the 
initial and flnal scale factor at the beginning and after inflation are related by Eflnai = W^^R iaj. It can be shown from 
equation (|37|) that the physical radius after inflation is given by ransi = lO^'^rinitiai. Thus we can see that upon substitution into 
equation (I76p . the 10*'^ factors in the numerator and denominator will cancel. Thus, it appears that Bt is also approximately 
constant during inflation, such that Bt will have approximately the same value before and after inflation. Therefore, after 
inflation one can continue to use equation (|74[) or (|76[) to express the cosmic primordial magnetic field. Then, after infiation, 
at the Hubble radius ~ 10"^'' cm, for t = 10"^" s, we find that ~ 2 x 10^^ gauss. 

Using equation (|74p or equation (|76|l we now calculate the cosmic magnetic field at recombination/decoupling and at the 
present epoch (-Bo). We will assume as suggested by particle physics that recombination/decoupling occurs ~ 350,000 years 
after the Big Bang (Liddle 2003). We can find values for this primordial frozen-in intergalactic (as related to spaces between 
galaxies) universal magnetic field Bt, as measured by an arbitrary conioving observer, at the present-day, normalized (see 
equations 1371 and I43p . spacetime coordinate distance rit = to) ~ 17 kpc [or z{t = to) ~ 4.6 x 10~^]. Recall that this observer 
is located at the centre of the Godel-Obukhov metric (equation [7]), which could be centred on an arbitrary comoving galaxy. 
Figure 6(c), shows how this cosmic Bt has evolved over time, at the above coordinate distance, from the spacetime separation 
r(t = 138 yr) to the present r(to). From this figure we see that Bt = -Brec ~ 4 x 10"^ gauss and Bt — Bo ^ 3 x 10~® gauss 
at recombination and at the present epoch, respectively. Importantly, the value for _Bo ~ 3 x 10~^ gauss is consistent 
with the upper limit constraint placed on the present strength of any primordial homogeneous magnetic field, which is 
-Bq ^ 4x 10~^ gauss, for fio = 1 and Ho = 71 kms~^ Mpc~^ (Barrow, Ferreira, & Silk 1997). The Cosmic Background Explorer 
(COBE) measurements provide this constraint: set by how the amplitude of the magnetic field is related to amplitude of the 
microwave background anisotropies on large scales. Pasquale, Scott, & Olinto (1999) study the effect of inhomogeneity on the 
Faraday rotation of light from distant quasi-stellar objects to find a consistent limit of _Bo ^ 10~^ gauss. 

The above present model calculated values of the cosmic magnetic field appear to be consistent with those that would 
allow the Universe to evolve into its present state, particularly like the cosmological model of Zeldovich (see Asseo & Sol 
1987, and references therein). Such cosmological models depend on the choice of the metric, of the equation of state, and 
of the cosmological constant Ao, where most often taken to be Ao = 0, until recently with the advent of dark energy. Such 
models are characterized by the expansion factors, the density of constituents, and the magnetic field, which evolve according 
to the Einstein-Maxwell field equations. The magnetic field energy is assumed to be modified only through the expansion 
process, and not through an exchange of energy with other constituents: matter or radiation, which means that the magnetic 
fiux is conserved and a uniform magnetic field evolves like 1/RxRy (for a magnetic field in the e^-direction), where Rx and 
Ry are scale factors. This results from the possibility for a uniform magnetic field to exist even in the absence of an electric 
current (V X B = 0). This can be interpreted as the evolution of the magnetic field compatible with the Einstein equations, 
as characteristic of a magnetic field geometrically frozen-in independently of the conductivity of the matter. 

Notice that the present epoch magnetic field of equation H74|) or H76p . goes to zero at the arbitrary comoving observer, just 
as does the gravitational acceleration (gGE)r of equation (|49|) . because r goes to zero. However, as we saw earlier in section[4l 
this is not the case for the GM acceleration {gGM)r of equation (|48p . which goes not go to zero, because the angular velocity 
iJiot of equation (|54p goes not go to zero. The former statement appears to be consistent with the relationship pointed out 
above between mass density and magnetic field through the spin density On the other hand, the latter statement is consistent 
with LJrot being related to both: {gGM)r and Bt, with Bt being related directly to the spin density (equation I73|) as well. 

5.6 The Godel-Obukhov Spacetime Equation of Motion Extension 

Now, if we assume that the theoretical model presented above is valid to describe processes in the Universe concerning the 
so-called dark energy, then the equation of motion of the scale factor of equation © can now be expressed as 
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in units of s~^, where, the GM acceleration of equation (|70|) . due to frame dragging, has been substituted in for the torsion 
term, which is the same as substituting equation ([Tip into equation ((5]). Note, 6, originally found in the metric of equation 
has been set equal to -j, as done in these model calculations (see equation I48|) . Equations (lllf) and (|28|) or (|48|) suggest 
that S = -J is the only well-behaved solution, i.e., for example, for z = 0.5, as measured by a present-day comoving observer, 
= gives the minimum value for ((?GM)r of equation (|48p . with this value increasing quickly to infinity at S ~ ±-|. It 
appears that S = f naturally sets up the geometry (in a rotating universe) for a local arbitrary comoving observer surrounded 
by concentric spacetime axisymmetrical hypersurfaces that look the same in radial directions along the line-of-sight of this 
observer, as discussed in Section [5.51 This is consistent with the quasi-Cosmological Principle of homogeneity and isotropy for 
such an observer who measures distances governed by, and from the centre of, the metric of equation Q. (Perhaps the use 
of cylindrical spacetime coordinates would clarify why 9 = ^ is a 'special' value.) Again, model parameters for k and a are 
defined in Section [3.51 B, as related to the spin density, is given above (equation [74] or [76} ; and m is given by equation ([4]). 
Note, the GM-inertial frame dragging term of equation (I77|) . believed, as suggested from these present model calculations, to 
be responsible for the recent accelerated cosmic expansion, and identified to be related directly to torsion of spacetime, is the 
only term dependent on the coordinate or spacetime separation r. This dependence on r is consistent with the accelerated 
cosmic expansion being observed to be a function of r (or z) as measured by a present-day arbitrary comoving observer, i.e., 
consistent with the apparently coming into being of the recent accelerated cosmic expansion at z « 0.5. As described by the 
GM acceleration of equation (I48p . this GM-inertial frame dragging (torsion) term decreases over time from large z to smaller 
z, but at some point in its later stage of evolution, it starts to increase with decreasing z, as can be seen in Fig. 2, and 
discussed in details in Sections 15.11 and 15.21 Thus, equation (|77|l appears to adequately explain the overall expansion of the 
Universe, wherein we consider cosmic rotation to be that of spacetime matter and is an inherent property of the geometrical 
and dynamical makeup of the Universe. 

Importantly, we see that the pressure p need not be negative to explain the recently observed cosmic acceleration, and 
Einstein's cosmological constant can remain his 'greatest blunder'; compare equation ((32}. We have found instead that this 
acceleration can be explained by considering the effects of the GM field and inertial spacetime frame dragging in a rotating and 
inertially expanding universe. So, in entirety, it appears that equation ((5}, equation ((77}, a possible extension of equation ((5}, 
and the metric of equation ((7} constitute a complete description of the spacetime dynamics of our expanding and apparently 
rotating Universe, past and present. 



6 CONCLUSIONS 

With the recent discovery of so-called dark energy (appearing to comprise ~ 70 per cent of the Universe), it seems we know 
little about the Universe we live in, save the mass-energy we can see and somewhat the effects of gravity this mass-energy 
displays. Now our understanding is limited even more in addition to our already lack of knowledge of what composes so-called 
dark matter (appearing to comprise ~ 30 per cent of the Universe). The above percentages are based on observations and 
the standard FLRW cosmology; and, importantly, they suggest that we should perhaps seek further an understanding of the 
underlining force: gravity. Such understanding could possible solve the problem of our recently, and somewhat embarrassingly, 
increeised lack of understanding. In this paper, I have adhered to the above suggestion by seeking an understanding of dark 
energy by considering it to be a component of gravity, a by product, arising in a general relativistically expanding and rotating 
cosmological spacetime. 

In this manuscript is presented a general relativistic model to describe the dynamical evolution of the Universe. This 
model appears to answered the question, 'Could dark energy be a manifestation gravity?' and the answer given is yes. In this 
model, the recently observed cosmic acceleration can be explained by considering the effects of the GM field due to inertial 
spacetime frame dragging in an expanding and rotating universe. These model calculations seem to show that application of 
the Godel-Obukhov spacetime metric of equation ((2}, or equation ((7}, in an Einstein-Cartan general relativistic spacetime, 
leading to the derivation of equation ((5}, the equation of motion of the scale factor R, and, perhaps, subsequently in the 
form of equation (I77|l , yield a complete description of how our Universe has evolved over time. Nonvanishing torsion of the 
Riemann-Cartan spacetime (Mao et al. 2007) and the GM field of inertial frame dragging appear to be one in the same in 
this description. Importantly, we see that the pressure p need not be negative, nor do we have to resurrect the cosmological 
constant A to explain so-called dark energy, as proposed by previous models constrained by the standard FLRW cosmology. 

Not only does the model presented here appear to describe the dynamical evolution of the Universe, but it appears to 
show a surprising relationship between the cosmic magnetic field and the GM field, through the spin density, which could 
possibly be a link between electromagnetism and gravity. This apparent relationship needs to be investigated further, not only 
on the astronomical level but the atomic level as well. However, this relationship should not be of a total surprise, since we 
have already made the analogy of the magnetic field of a rotating charge and the GM field of a rotating mass, based on their 
similarities, and, particularly, how the Lorentz force on a moving charged particle in a magnetic force field is analogous to the 
Coriolis-like force on a moving particle in a GM force field (see e.g., Braginsky, Caves, Thorne 1977; Thorne et al. 1986). 
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As a future calculation, one might use the results of the cosmological model presented in this manuscript to evaluate 
the Friedmann-like equation, describing the evolution of the scale factor, derived from the gravitational field equations by 
Obukhov (2000), used in the derivation of equations ([5]), (|7ip . and (|77|l . It appears that equation (5.42) of Obukhov (2000) 
can be used, applying the Equivalence Principle, to see what fraction of the critical density the GM acceleration, due to frame 
dragging (or torsion), contributes to making the observed density parameter SI ~ 1. This might also shed light on the true 
nature of dark matter. 

Finally, whether all the concepts presented in this manuscript are fully valid, as they appear, or not, particularly those 
leading directly to equation (|77|l . the fact remains that the Godel-Obukhov cosmology of equation ((Sj exhibits cosmic accel- 
eration in the equation of motion of the scale factor, which could possibly explain the recently observed cosmic accelerated 
expansion. If the extended form (equation 1771) is indeed valid it shows that this cosmic acceleration became important in 
recent times. 
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